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The Harman method for measuring the thermal conductivity of a sample using the Peltier effect, 
may also be used to determine the dimensionless figure of merit from just two electrical resistance 
measurements. We consider a modified version of the Harman method where the current contacts 
are much smaller than the contact faces of the sample. We calculate the voltage and temperature 
distributions in a rectangular sample of a material having anisotropy in all of its transport coef- 
ficients. The thermoelectric anisotropy has important consequences in the form of thermoelectric 
eddy currents and the Bridgman effect. We prove that in the limit of a very thin sample of arbitrary 
shape, there exist van der Pauw formulae relating particular linear combinations of the potential and 
temperature differences between points on the edges of the sample. We show that the Harman figure 
of merit can be radically different from the intrinsic figures of merit of the material, and can often 
be substantially enhanced. By defining an effective figure of merit in terms of the rate of entropy 
production, we show that the increase in the Harman figure of merit does indicate an improvement 
in the thermoelectric performance of an anisotropic sample having small current contacts. However, 
we also prove that in the case of a material with tetragonal symmetry, the effective figure of merit 
is always bounded from above by the largest intrinsic figure of merit of the material. 
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I. INTRODUCTION 

It is widely accepted that improving the efficiency of 
thermoelectric devices is one of the key challenges cur- 
rently facing both technologists and physicists alike. A 
thermoelectric device either converts heat into electrical 
power through the Seebeck effect, or employs the Peltier 
effect to transport heat. The efficiency of such devices 
depends primarily on the transport coefficients of the 
constituent thermoelectric materials. The ideal thermo- 
electric material would have i) a large thermopower S 
to give maximum transfer of heat into electrical power or 
vice versa; ii) a large electrical conductivity a to minimise 
the wasted power associated with Joule heating; and iii) a 
low thermal conductivity k to minimise thermal shorting. 
The thermal conductivity will, in general, include both 
contributions from the lattice, and from the charge carri- 
ers themselves. The actual efficiency of a thermoelectric 
device will generally depend in a very complicated way on 
the transport coefficients of its composite materials [![. 
However, a simple result is obtained in the case of a single 
thermoelectric material operating as a generator between 
two temperatures, and whose transport coefficients are 
taken to be temperature independent. The maximum 
efficiency of such a device turns out to depend on the 
transport coefficients only through a quantity known as 
the figure of merit (FOM), Z — S 2 a/n. The larger the 
value of Z, the higher the efficiency. The focus is there- 
fore on discovering or designing materials with large fig- 
ures of merit. A dimensionless FOM, z = ZT in excess 
of unity is needed to make thermoelectric devices com- 
petitive in terms of efficiency ■ 

The ability to accurately determine the FOM is of 
great importance. One way to determine the FOM is 
simply to make separate measurements of all three trans- 
port coefficients individually. However, there is a com- 



monly employed technique which allows a determination 
of the dimensionless FOM from just two electrical re- 
sistance measurements. This is known as the Harman 
method [H, and was originally designed as a means of 
measuring the thermal conductivity using the Peltier ef- 
fect. An electrical current passed through the sample re- 
sults in a temperature difference across the sample due to 
the Peltier effects occurring at each contact. The steady 
state temperature difference results from the balance be- 
tween the Peltier effects and the heat flow by thermal 
conduction, which allows the thermal conductivity to be 
determined. It turns out that the electrical resistances of 
the sample due to alternating and direct current in the 
Harman set-up allow the dimensionless FOM of the sam- 
ple to be determined from z = Rdc/ R&c — 1- This result 
depends on a large number of assumptions, mainly relat- 
ing to the various heat losses and heating effects which 
can occur in the sample. However, there is an additional 
assumption, which is that the flows of heat and electricity 
are uniform through the sample. It is important to ask 
what would happen if this assumption were not fulfilled. 
Such a situation might occur if the contacts between the 
sample and the current leads are smaller than the cross- 
sectional area of the sample itself. 

We will show later in this work that the effects of the 
small current contact do not change the FOM measured 
with the Harman method provided that the transport 
coefficients of the sample are isotropic. However, it is 
extremely unlikely that this would occur in practice. In- 
deed, many of the materials which are currently under in- 
vestigation as potential thermoelectrics have anisotropic 
transport properties. These include the layered cobalt 
oxides [4], and semiconductor superlattices [E S 0j- 
Even traditional materials such as BiaTea have uniax- 
ially anisotropic transport properties [8j. In particular, 
anisotropy in the thermopower can lead to new and inter- 
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esting effects not seen in isotropic systems. First, there 
is an additional internal heating effect, the Bridgman ef- 
fect 0, [l(| , which occurs if the current distribution be- 
comes non-uniform. Importantly, this effect is linear in 
the electrical current, and hence cannot be neglected as 
a higher-order process. Second, the steady state with 
an applied temperature distribution generally contains 
circulating electrical currents, known as thermoelectric 
eddy currents [ll|. These lead to a non-trivial relation- 
ship between the voltage and temperature distributions 
in the sample. 

The theory of electrical resistance measurements with 
small current and voltage contacts is extremely well- 
established in the literature. The classic work of van der 
Pauw [l2j considered electrical resistance measurements 
in thin samples of arbitrary shape. It was shown that two 
particular four-terminal resistance measurements on an 
isotropic sample obey a universal relationship, allowing 
the electrical resistivity to be determined. Logan et al 
[l3| considered the problem of determining the voltage 
distribution in a rectangular sample of isotropic resis- 
tivity with current contacts on the corners. Their solu- 
tion was subsequently used by Montgomery [l4[ to show 
how all three components of the resistivity tensor in an 
anisotropic sample can be determined from appropriate 
resistance measurements. To our knowledge, the equiva- 
lent calculations have not been extended to thermal and 
thermoelectric measurements. 

A motivation for this work is recent experimental work 
[l5l HU where the Harman method is used to measure 
the FOM of various anisotropic materials. The contact 
configurations used in the experiments are very similar 
to those used in Montgomery's methods of measuring 
the electrical resistivity, with small current and voltage 
contacts on the surface of the sample. In some configu- 
rations, the Harman FOM (i.e. the FOM measured us- 
ing the Harman method) is in good agreement with the 
known intrinsic FOM of the material. However, in mea- 
surements on (Bi .8Pbo.2)2Sr 2 Co20y (BiPbSrCoO) sin- 
gle crystals in the out-of-plane direction, the Harman 
FOM is found to be substantially enhanced above the 
intrinsic FOM. A key feature of this material is that 
while the thermal conductivity is relatively isotropic, the 
electrical conductivity exhibits an anisotropy of around 
10 4 between in and out-of-plane directions. We will dis- 
cuss the results of these experiments in some detail later 
in the paper. We also note that other experiments on 
(Bi, Sb)2Te3 [l?]] have also revealed an intriguing depen- 
dence of the Harman FOM on the sample length, cross- 
sectional area and voltage contact location. 

The principle aim of this paper is to calculate the 
voltage and temperature distributions in a sample with 
small current contacts, and having arbitrary degrees of 
anisotropy in its transport coefficients. We will then in- 
vestigate how this alters the FOM measured in the Har- 
man method. We will show that the Harman FOM no- 
longer takes a single value for the sample, and is in gen- 
eral a function of the points between the resistances are 



measured. In some configurations the Harman FOM is 
substantially enhanced above the intrinsic FOM of the 
material in the measurement direction. By defining an 
effective FOM based on the rates of entropy production 
in the Harman method, we show that the increased Har- 
man FOM does indicate an improvement in the thermo- 
electric performance of the system. We identify two im- 
portant mechanisms responsible for this increase. First, 
in a system whose electrical and thermal anisotropics are 
different, heat and electricity can flow along very differ- 
ent paths. Second, thermoelectric eddy currents and the 
Bridgman effect lead to a mixing of the thermopowers 
in the different coordinate directions, and can lead to an 
increase in the size of the thermoelectric effects. We will 
also prove, however, that the effective FOM is always 
bounded from above by the largest intrinsic FOM in an 
anisotropic material. 

The rest of this paper is organised as follows. In sec- 
tion [Hi we set-out the theory needed to describe trans- 
port in anisotropic systems. In particular, we will discuss 
the origins of the new thermoelectric phenomena, namely 
eddy currents and the Bridgman effect. In section UTTI wc 
describe the Harman method in some detail, discuss the 
various assumptions made, and show how some of the 
errors which occur can be reduced. In section IIVI we 
present the model, and obtain the most general solution. 
The solution can be simplified once the current contacts 
locations are fixed, and we derive the simplified solutions 
for two different configurations in appendices and [B] 
The Harman FOM is calculated in section IVT1 for several 
different cases of anisotropy and sample geometry. In ad- 
dition, we consider the limit of a very thin sample, and 
show that equivalents of van der Pauw's relation hold in 
the Harman method. In section I Villi we apply the re- 
sults to the recent experiments of Kobayashi et al and 
Tamura et al. Finally, in sections ILXI and [X] we inves- 
tigate whether the increased Harman FOM indicates an 
improvement in thermoelectric performance. This leads 
us to a new definition of the FOM in terms of the rates 
of entropy production in the Harman method. 



II. THEORY OF ANISOTROPIC TRANSPORT 

In an anisotropic system, the zth components of the 
electrical and heat current densities, ji and jq, are 
given respectively by the linear response equations (xi — 
{x,V,z}) 



]i — &ik- 
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where the summation convention has been invoked. Here, 
T = T(x,y,z) is the absolute temperature, and x = 
x(x,y,z) is the electrochemical potential. This is de- 
fined as x — A*/ e ~~ 0j where /1 is the chemical potential, 
(j> is the electrical potential, and e is the magnitude of 
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the electronic charge. The Kelvin/Onsager relations im- 
ply that the Peltier and Seebeck coefficients are related 
by Uik = TSki- In this work, all of the transport coef- 
ficients are taken to be diagonal tensors. Effects associ- 
ated with the off-diagonal components [H, Oil will not 
be considered here. 

Let us now consider the Seebeck effect in an anisotropic 
system. In the conventional picture of a homogeneous, 
isotropic system, a temperature gradient applied to an 
isolated sample causes an electrical current to flow, until 
the build-up of charge creates a voltage gradient which 
counteracts it. This leads to a steady state in which no 
net electrical current flows. Setting ji = and Sik — 
SSik in equation ([I]), and solving for \ gives the familiar 
result that the thermoelectric voltage is given by (up to 
an arbitrary constant of integration) xth — ST in the 
steady state. 

When the thermopower is anisotropic, the steady state 
does not have ji — in general. To see this we first 
rearrange equation (JTJ) for dx/dxi 



dxi 



Pikjk + Sik 



dT 

dx k ' 



(3) 



where p = a^ 1 is the resistivity tensor. If ji = is 
a solution of this equation, then we should be able to 
integrate to find x m terms of T. This is only possible if 
Vx is a conservative vector field, and hence V x (Vx) = 0. 
Setting ji — and requiring that the curl of the right- 
hand side of equation ([3]) equals zero implies that 



(S X Sy) 



d 2 T 
dxdy 



q2 T q2 T 
0=(S x -S z )i-^- = (S v -S z )- 



dxdi 



dydz ' 



where S x = S xx etc. Hence unless the temperature dis- 
tribution depends on only one coordinate, we see that 
the steady state where all components of the current 
density vanish cannot exist when the thermopower is 
anisotropic. This implies that in the steady state there 
are electrical currents flowing in the bulk, providing an 
Ohmic component to the voltage which ensures that 
V x (Vx) = 0. These currents, often referred to as 
thermoelectric eddy currents, were first predicted and 



observed by Samoilovich and Korenblit [11|. We note 
that their existence has the following important conse- 
quence; the relationship between the thermoelectric volt- 
age and temperature distributions does not generally obey 
Xth = ST when the thermopower is anisotropic, even if 
S itself is independent of temperature. 

We now turn to the internal heating effects which can 
occur in an anisotropic system. These effects may be de- 
rived by applying the principles of conservation of charge 
and conservation of energy |10j , which can be expressed 
in terms of continuity equations dp/dt + V.j = and 
du/dt + V.jy — 0. Here, p and u are the charge and 
energy densities respectively. The quantity ju is the en- 
ergy current density, related to the heat current density 
via ju — jq — xj- In a steady state, the continuity equa- 



tions become 



dj 



dji 
dxi 



0. 



(4) 



m _ df(j,g ik _ x s ik )j k - Kik ^L\ = o. (5) 
oxi axi V oxk ' 

Expanding out the right-hand side of the second of these 
equations and making use of equation |T]) leads to the 
following equation: 



dT 



d ( dT\ 

-7T— I Kik q = Plkjljk 

axi V oxk ' 
rr . fdS lk \ f dS ik \. dT dj k 



The left-hand side is the divergence of the heat current 
which would occur in the absence of an electrical cur- 
rent. The terms on the right-hand side are source terms 
which represent the internal heating effects. The first 
term is recognised as the Joule heating, while the remain- 
ing terms contain the thermoelectric effects. The first two 
terms on the second line are the (internal) Peltier effect 
and the Thomson effect, both of which are familiar from 
isotropic systems. These terms originate from the fact 
that if the thermopower varies with position, then as an 
electrical current flows through the sample the amount of 
heat it carries also changes. This variation can come from 
an intrinsic inhomogeneity of the sample (Peltier effect), 
or because the thermopower is a function of tempera- 
ture which itself varies throughout the sample (Thomson 
effect). The most common manifestation of the Peltier 
effect occurs at the junction between two different mate- 
rials, where the rate of heating at the junction is given 
by the difference of the Peltier coefficients in each ma- 
terial multiplied by the current, Q = (H2 — If 
the thermopower is anisotropic, then we must also take 
into account the direction in which the current enters the 
sample, since this will determine which component of the 
Peltier tensor appears in this expression. 

The third term on the second line represents an effect 
originally proposed by Bridgman [9]. Unlike the Peltier 
and Thomson effects, the Bridgman effect only occurs 
when the thermopower is anisotropic. For isotropic ther- 
mopower Sik = SSik, and the Bridgman term reduces 
to —TSdji/dxi, which vanishes in the steady state. To 
understand the physical origin of this term, consider a 
current flowing in the x-direction (say) in an anisotropic 
material. It carries an amount of heat proportional to 
IIj;. If the current were to change direction, say into 
the y-direction, it would now carry a different amount 
of heat, now proportional to n y . The Bridgman effect 
therefore captures the heating and cooling effects asso- 
ciated with changes in the direction of the current in a 
system with an anisotropic thermopower. 

Having outlined the basic theory we shall require in 
this paper, we now turn to a discussion of the conven- 
tional Harman method. 
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III. THE HARMAN METHOD 

In the 1950's, Harman [|[ proposed a novel technique 
for measuring the thermal conductivity of a material, 
which utilises thermoelectric effects. In the conventional 
set-up, current leads are attached to the sample in such 
a way that the electrical current flows in the direction in 
which the thermal conductivity is to be measured, the 
^-direction say. The leads are then used to suspend the 
sample in a vacuum. As the current enters sample, the 
Peltier effect causes a heating/cooling of the contact area, 
with an equal cooling/heating as it leaves through the 
other contact. The rate of heating/cooling is given by 
(II; — Yl z )I, where 11/ and H z are the Peltier coefficients 
of the leads and the sample (in the z-direction) . The tem- 
perature gradient created in the sample results in a back- 
flow of heat by thermal conduction, which eventually bal- 
ances the Peltier effects leading to a steady state. The 
simplest treatment of this problem makes the following 
assumptions; (i) All internal heating effects (as described 
in section[TT| in the sample may be neglected; (ii) All heat 
losses from the sample may also be neglected; (iii) The 
electrical current is injected/removed uniformly across 
the contact faces of the sample. Hence, the Peltier effects 
cause the contact faces to be heated/cooled to uniform 
temperatures Th/T c , and the flows of heat and electricity 
may be treated as being effectively one-dimensional. 

Under these assumptions, it is straightforward to solve 
equation ([6]) for the temperature distribution in the sam- 
ple, 



T(x,y,z) = T- 



n 7 ab 



n z )Icrz 



(7) 



where the contact faces of the sample lie at z = 0, c, and 
have an area a x b. The ambient temperature is denoted 
by T. It follows that the temperature difference between 
any two points (x, y, 0), (x',y',c) is 



AT = T(x',y',c) 



-T(x,y,0) 
i/c _ (Ih 
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where K z q = n z ab/c is the thermal conductance of the 
sample in the z-direction. If II; and Tl z are already 
known, then measuring AT determines k z . 

In addition to measuring the thermal conductivity, the 
experimental setup of the Harman method may also be 
used to measure the FOM of the sample directly The 
voltage distribution in the sample with a direct current 
flowing can be found from equation Q: 



X(x,y,z) =x- 



Ic (z 1\ S z (Ui - n z )Ic (z I 



ab\c 2 



K z ab 



(9) 

The potential difference measured across the sample us- 
ing leads with a thermopower Si is given by 



V dc = IR zQ - [Si - S Z )AT. 



(10) 



In a second measurement, the potential difference is again 
measured across the sample, but now with an alternating 
current of amplitude / flowing. Provided that the period 
of the current is much shorter than the characteristic time 
for the temperature difference to develop, there will be 
no contribution from the thermoelectric term, and hence 

Kc = IR z0 , (11) 

where R z q = c/a z ab is the electrical resistance of the 
sample in the z-direction. Defining a quantity R dc = 
Vdc 1 1 and using equation ([5]), the potential difference 
can be written as 



IRdc = I Rzo + (Si — S z ) 



(Tli-Tlz)I 

K z0 



(12) 



Using U l -U z = T(Si -S z ), the FOM as determined by 
the Harman method is then given by 
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(13) 
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Under 'normal' circumstances the product R z qK z q — 
p z K z , since the area and length factors cancel out. The 
FOM measured by the Harman method is a relative FOM 
between the sample and the measuring leads, since the 
thermopower appearing in (fT3| is Si — S z . The ther- 
mopower of the leads is typically much smaller than that 
of the sample, and so (Si — S z ) 2 S z . Hence this Har- 
man FOM is approximately equal to the intrinsic FOM 
of the sample (in the z-direction) . 

Let us now consider the validity of the assumptions 
(i)-(iii), which were used to derive the results above. As- 
sumption (i) is extremely crude, and requires further 
investigation. By restricting attention to homogeneous 
samples, we can eliminate the internal Peltier effects, 
while the Thomson effect is small and can usually be ne- 
glected. The Bridgman effect is also eliminated provided 
that assumption (iii) holds, since the current density is 
uniform throughout the sample. Note that the Id na- 
ture of the temperature distribution also guarantees that 
there are no thermoelectric eddy currents in the sample, 
and that the thermoelectric voltage and temperature dis- 
tribution are related by Xth = ST. A related assumption 
is that the temperature difference generated is sufficiently 
small that linear response equations may be applied; this 
can be achieved by carefully choosing the current. 

The effects of Joule heating are relatively easy to in- 
clude within the one-dimensional model, since the rate of 
heat production is constant throughout the system due 
to the constant current density. A finite thermal conduc- 
tance in the leads needs to be included to allow the Joule 
heat to be removed. In addition, we include a contact re- 
sistance rt and for each current contact, which may 
be different. Again solving equation ((6|), we find that the 
temperature difference between the contact faces of the 
sample is now given by 



AT = T(x',y',c)-T(x,y,0) 
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where Ki is the thermal conductance of a single lead. 
The result is completely independent of the bulk resis- 
tance, and hence we conclude that provided the tem- 
perature is measured between the extreme edges of the 
sample, the bulk Joule heating has no direct effect on 
AT. The effect of conductive heat loss along the leads 
reduces the temperature difference across the sample, as 
one might naively expect. The current leads will often 
have a much lower thermal conductance than the sam- 
ple, and this term can be ignored. If the contact resis- 
tances are not equal, then a convenient way to eliminate 
their effect is to make measurements of AT for currents 
passed in both directions through the sample, and cal- 
culate AT a = (AT (I) - AT(-7))/2. This removes the 
contact resistance contribution, since it is proportional 
to I 2 ||. 

We have already dealt with one aspect of assumption 
(ii) above, namely the heat loss due to thermal conduc- 
tion along the current leads. Since the sample is hung 
in vacuum, convective heat losses are eliminated, leaving 
radiation as the only other possible mechanism. Heat 
losses due to radiation effectively leads to an additional 
flow of heat through the sample, and makes the thermal 
conductivity appear too large. Appropriately choosing 
the aspect ratio of sample can help to reduce this error 




The final assumption (iii) that the flows of heat 
and electricity may be treated as being effectively one- 
dimensional, does not seem to have been considered pre- 
viously. The one dimensionality relies the contacts being 
made in such a way that the temperature and voltage 
are uniform across the contact faces of the sample. In 
Harman's original work it was noted that this could be 
achieved by applying a thin layer of solder over the con- 
tact faces of the sample. The high electrical and ther- 
mal conductivity of the solder ensures assumption (iii) is 
valid. However, one could ask what would happen if this 
condition were not satisfied, either through error (i.e. the 
solder does not completely cover the surface of the sam- 
ple, or does not make a good contact in some places), or 
by deliberately using small contacts. It is reasonable to 
assume that this will lead to more complicated distribu- 
tions of temperature and voltage in the sample. If the 
thermopower is anisotropic, we have already seen that 
new effects can occur when the temperature and volt- 
age distributions become non-trivial, both in the form of 
eddy currents and the Bridgman effect. Even if the ther- 
mopower is isotropic, we will see that electrical and ther- 
mal anisotropy can also have interesting consequences. 
A further degree of freedom is the geometry of the sam- 
ple itself, and the particular arrangement of current and 
voltage contacts on its surface. In the next section, we 
introduce the model which will be used to investigate the 
effects of small contacts and anisotropy in the Harman 
method. 
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FIG. 1: Schematic representation of the sample geometry. 
All faces are electrically and thermally insulated. An electri- 
cal current / is passes into the sample through the contact at 
(asi, 2/1, zi), and leaves the sample through an identical con- 
tact at (x2, J/2, Z<i). (a) Point current contacts; (b) Current 
contacts with a finite area d 2 . 



IV. MODEL 

We seek a minimal model which is capable of exploring 
the effects of small current contacts and transport coef- 
ficient anisotropy on the temperature and voltage distri- 
butions created in the Harman method. We begin with 
a description of the sample configuration. The sample is 
chosen to be a cuboid with dimensions a x b x c in the 
(x, y, z) directions respectively. We choose a coordinate 
system which coincides with one corner of the sample, 
as illustrated in figure [1] The sample material is taken 
to be homogeneous, but may have arbitrary degrees of 
anisotropy in all of its transport coefficients. 

The sample is connected to two current leads, which 
have a thermopower 5;, and which we assume to have 
negligible thermal conductance. Hence, there is no heat 
loss from the sample along the current leads. The leads 
connect to the sample through square contacts of dimen- 
sions d x d. Each contact may be placed on any face 
of the sample. The boundary conditions at the current 
contacts need careful consideration. In reality, the con- 
tact is likely to be made from a blob of solder, which 
will have a very high electrical and thermal conductivity. 
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The region under the contact will therefore be approxi- 
mately equipotential and isothermal. The electrical and 
heat current density will vary under the contact to en- 
sure that these conditions are satisfied. Away from the 
current contacts, all surfaces of the sample are taken to 
be electrically and thermally insulated. The disadvan- 
tage of these boundary conditions is that they make the 
problem very difficult to handle analytically. The reason 
is that the boundary conditions on the contact faces are 
not uniform, since over one part of the surface the func- 
tion itself is specified, while over the rest the derivative 
is specified. To avoid this complication, we adopt a dif- 
ferent boundary condition, namely that the normal com- 
ponents of the electrical and heat current density under 
the contacts are uniform. This means that the problem 
has Neumann boundary conditions on all surfaces, and 
can be solved by separation of variables. In the limit 
of point contacts, the two choices of boundary condition 
are in some sense equivalent, although the voltage and 
temperature are infinite at the point contacts. For small 
contacts (i.e. much smaller than the cross-sectional area 
of the sample) , the choice of boundary conditions is not 
expected cause large changes in the solution in parts of 
the sample away from the contact. Similarly for large 
contacts (approximately the same area as the sample) 
the variation of the temperature and voltage in the con- 
tact region is expected to be small. Finally, we will show 
later although the variation of the temperature and volt- 
age under the contacts is unphysical, the average tem- 
perature and voltage underneath the contacts does have 
a physical significance. 

Since we require a minimal model only, we will not 
attempt to solve the formidable non- linear problem posed 
by equations ^ and ©■ Instead, take the following 
approach to linearise these equations. We first write the 
temperature and electrical current density in terms of 
symmetric and anti-symmetric parts with respect to the 
injected current I; j(J) = j s (I) + j a (I), T(I) = T S (I) + 
T a (I), with j s (-7) =F(I), 3 a (-I) = -j a C0 etc. Now 
consider expanding j and T in powers of /. Clearly, we 
have j a (I)_ = 0(1), j s (I) = 0(I 2 ), T a (I) = <D{I) and 
T S (I) = T + 0(I 2 ). Substituting these expressions into 
equations (UJ and ([6]), and retaining only terms of O(I), 
we find 



with 



= 0, 



d f dT a 

oxi V axk 



-TS. 



ik 



dxi '' 



(15) 
(16) 



d X a q dT a 

OX k OX k 
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All transport coefficients are to be evaluated at the ambi- 
ent temperature T . Note that the internal Peltier effect 
is absent, since we have restricted attention to homoge- 
neous samples. The Thomson effect and the Joule heat- 
ing also do not appear, since they are both 0(I 2 ) effects. 
Further, they would only alter the anti-symmetric parts 
of x an d T at 0(I 3 ). Hence, by working with T a and x" 'j 
our approximation holds to 0(I 3 ). In the remainder of 
this paper we will drop the 'a' superscript, but it should 
be remembered that when comparing with experiments 
the anti-symmetric parts of x and T should be taken. As 
discussed in section ITTA this procedure is usually adopted 
in experiments anyway in order to remove extraneous ef- 
fects of Joule heating associated with contact resistance 
or misplacement. 



V. SOLUTION 
A. General Solution and the Reciprocity Theorem 

We will now find the most general solution to equa- 
tions ([4]) and ([6]). In calculating the voltage distribution 
in an isotropic sample with point electrical contacts on 
the corners, Logan et al [l3[ used an analogy with elec- 
trostatics to obtain the solution. Each point contact acts 
as a source of current, much like a point charge acts as 
a source of electrical flux. An infinite number of im- 
age charges must be introduced to enforce the boundary 
conditions on the surfaces of the sample. The voltage 
distribution in the sample is then given by the sum of 
the potentials due to each image charge. 

We will adopt a similar method to solve equations (fT5|) 
and (fl6|) . To illustrate the method, we consider two 
point current contacts located at positions (xi,yi,zi) 
and (x2, ?/2, z 2 ). I n practice, these points must lie on the 
surface of the sample, but for argument sake we consider 
them to lie in the interior of the sample. One contact is 
a source of current, while the other is a sink. Defining 
A(x, y,z) = 6(x- x 1 )8(y - y\)5(z - zi)-S(x- x 2 )5(y - 
2/2)^(^ — ^2), equations (fT"5|) and (|16p . now written-out in 
full, are modified to read 



d 2 x ^ d 2 X ^ d 2 X Q d 2 T 
ox z ay A oz z ox z 



d 2 T 



d 2 T 



(18) 



(k x +T(j 



Q 2^d 2 T 2 d 2 T 2 ,d 2 T - ( 5 2 y „ d 2 x q d 2 X 

•^J^-(- y +^^)^-(^+^^)^+^(-^^+^^^2+^^^2 



ILiIA(x,y,z). 

(19) 
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The right-hand side of equation (|18[) represents a point 
source of current at (x±,yx, Zi), and a point sink of cur- 
rent at (x2, yi,Z2)- The right-hand side of equation (fl9|) 
represents a point source of Peltier heat at {x\,y\,Z\), 
and a point sink of Peltier heat at (£2,2/2,22)- Rather 
than attempt to explicitly find the locations of the image 
charges, we adopt a different approach. We expand \ 
and T using a basis of functions whose derivatives vanish 
at the sample boundaries: 

X = /J dmnp cos fhx cos ny cos pz, 

mnp 

T = b mnp cos ma; cos hy cos pz, 

mnp 

with to = mix /a, h = nir/b and p — pix jc, and where 



the sums run over all integers, —00 < to, n,p < 00. It is 
clear from equations ^ and §2§ that by requiring that 
the derivatives of x an d T vanish independently on all 
boundaries, the sample is electrically and thermally in- 
sulated. Similarly, A(x,y,z) may also be written in this 
way, using the following representation of the delta func- 
tion [12|: 

S(x — x\) = — 'S~ y cosrhxi cosrhx. (20) 



We now substitute these representations into equations 
(fT8| and (fT9| , which reduce to two simultaneous algebraic 
equations for a mnp and b mnp . Omitting the algebraic 
details, we find for x and T 



1 



x(x,y,z) = 1^2 



cos UIX2 cos ny 2 cos pz 2 — cos mx\ cos ny\ cos pz\ 



abc 

mnp 



■ cos Tnx cos ny cos pz x 



{k x - VxSJJh - ^x))m 2 + (k v - <7yS y (Ui - U y ))h 2 + (k z - a z S z (Ui - U z ))p 2 

(a x m 2 + o y h 2 + a z p 2 )(K x m 2 + n y h 2 + K z p 2 ) + f(a x a y (S x - S y ) 2 rh 2 fi 2 + o x o z (S x - S z ) 2 m 2 p 2 + a y a z {S y - S z ) 2 h 2 p 2 ) 

(21) 

. . cos mx2 cos ny 2 cos pz2 — cos mxi cos ny 1 cos pzi _ 

l(x,y,z) = —l y cos vcix cos ny cos pz x 

' abc 

mnp 



a x (Ili - U x )m 2 + a y (Ui - U y )h 2 + cr z (n ; - U z )p 2 



(a x m 2 + Oyh 2 + a z p 2 )(n x m 2 + n y n 2 + n z p 2 ) + T(a x a y (S x - S y ) 2 m 2 n 2 + o x o z (S x - S z ) 2 rh 2 p 2 + a y a z [S y - S z ) 2 h 2 p 2 ) 

(22) 



The prime above the sum indicates that the term with 
in = n = p = is not included. Physically, this is be- 
cause there is no net current /heat being delivered to the 
sample. Although these expressions may seem somewhat 
cumbersome, the elegance of this representation is that 
it emphasises the symmetry of the solutions with respect 
to the three coordinate directions. Only the choice of the 
current contact positions breaks this symmetry. 

Another advantage of this representation is that it en- 
ables us to easily investigate a property of the solution 
known as reciprocity. Define the quantity Vabcd as 
the potential difference measured between points C = 
(x,y,z) and D = (x',y',z'), with a current Iab en- 
tering the sample at A and leaving at B. The reci- 
procity theorem states that the voltage is symmetric 
under interchange of voltage and current contacts, or 
Vabcd = Vcdab, provided that total current passed 
is the same, i.e. Iab = Icd- This theorem is known to 
hold for the ac voltage i. e. in the absence of thermoelec- 
tric effects. The form of equation (f2"Tj) indicates that the 
theorem continues to hold in the presence of (anisotropic) 
thermoelectric effects. Further, it is clear from equation 



(|22| that this reciprocity also holds for the temperature 
difference, 1ST abc d — ^Tcdab- That reciprocity con- 
tinues to hold in this problem is not surprising, given that 
the linear response condition required for its existence is 
satisfied here (23j . 

Once the particular faces on which the contacts are lo- 
cated have been chosen, it is possible to simplify the form 
of solution by performing one of sums appearing in (|2ip 
and (f2"2"]) . There are three distinct contact configurations 
on a cuboid. These have the contacts on the same face, 
on opposite faces, and on adjacent faces of the sample. In 
this work, we will primarily be interested in the first two 
of these configurations. The choice of particular faces for 
the calculations is of course arbitrary; here we will choose 
the faces to to those at z = 0, c. Although the solution 
with contacts on adjacent faces follows naturally from 
equations (|2"Tj) and (|2"2"]) with appropriate choices of con- 
tact location, it is not easy to simplify the solution into a 
compact form. Hence, we will not explicitly consider this 
configuration in this paper. The solutions for contacts on 
the same side and on opposite sides of the sample may be 
found in appendices [AlandlBl In the next section, we will 
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explore the general form of the solution by considering 
the voltage distribution due to an alternating current. 



B. The ac Voltage 

A quantity of particular interest for each contact con- 
figuration is the voltage distribution in the sample due 
to an alternating current. As discussed in section IIII1 
the Peltier effects at the current contacts are unable to 
change the temperature distribution in the sample if the 
current oscillates with sufficiently high frequency Since 
the temperature in the sample remains uniform (at lead- 
ing order in I), all thermoelectric effects vanish. The ac 
voltage distributions for each contact configuration are 
obtained from equations (|A8j) and (|B2|I by setting all 
thermoelectric coefficients equal to zero: 

Ic z 1 I ' 

mn 

(23) 

Xac = -TJ 22 ^mn(x, V, Z\ R m n) , (24) 

a z d l ^ 

ran 

where F^ nn = r 2 z m 2 + r 2 z h 2 . Here, Xac i s the ac volt- 
age distribution in a sample whose current contacts are 
on the same face at z = 0, while xlc 1S the ac voltage 
distribution in a sample whose current contacts are on 
opposite faces at z = and z — c. The current contacts 
each have a finite area d 2 . The solution for x! c divides 
into two parts. The first term corresponds to a uniform 
flow of electrical current through the sample. The sum 
represents the modification to this uniform flow due to 
the small current contacts. The ac voltage has a mini- 
mum over the contact where the current enters the sam- 
ple, and a maximum over the contact where the current 
leaves. The maximum (minimum) values of the voltage 
occur at the centers of the current contacts, and increase 
(decrease) as the contact size d tends to zero. The scale 
over which the voltage distribution is modified in the z- 
direction is controlled by the electrical anisotropies r xy 
and r yz . Larger anisotropy leads to the voltage distribu- 
tion being modified only very close to the contact faces. 

Very similar statements can be made for xtc- The main 
difference from the previous case is that there is no linear 
term in the solution. This is simply because there is no 
uniform flow of electricity in the ^-direction; the current 
enters through one contact on the face at z = 0, and 
leaves through another contact also at z = 0. 

VI. HARMAN FIGURE OF MERIT 

In section IIII1 we showed by direct calculation that 
the Harman FOM given by equation (fT"3"]) is equal to 
T(Si — S z ) 2 cr z / 'k z , which is approximately equal to the 



intrinsic FOM in the z-direction, z z . Here, we will sim- 
ply define the Harman FOM by the equation zh = 
(Vdc — V^c)/\4c and compute this quantity using the re- 
sults of the previous section. 

It is obvious that once the temperature and voltage 
distributions vary with position on the contact faces, the 
Harman FOM is no-longer a unique quantity. In general, 
it is a function of twelve coordinates; the two current 
contact locations, and the two voltage contact locations. 
It is also generally true that the magnitude of the Harman 
FOM will be equal to zero for certain configurations, and 
tend to infinity for others. This simply because for given 
contact contact locations it will generally be possible to 
find voltage contact locations for which either Vd c — V ac = 
with V„ i= 0, or V„ = with V dc - V ac ^ 0. 

Rather than state the results in full generality, it is 
instructive to start with the simplest case first. 



A. Large Contacts 

We start with the case where the sample has a square 
cross-section (a = b), and the contact size is made equal 
to the area of the sample, d = a. Clearly, this case is 
only applicable when the contacts are on opposite sides 
of the sample. Setting d = b = a, we find that all of 
the Fourier coefficients (|A2| apart from j/oo and 700 van- 
ish identically. This means that the sums appearing in 
equations (|B2|) are also identically zero. Only the lin- 
ear terms remain, the expressions for \° and T° being 
identical to © and (J7J) respectively, up to the constant 
mean temperature and voltage which we have neglected. 
Hence the temperature and voltage are both uniform 
over the contact faces. Following through the analysis 
presented in section IIII1 the Harman FOM reduces to 
zh = T(Si — S z ) 2 a z / k z , which is the expected value. 
Note that this result does not require any assumptions 
about the anisotropy; large current contacts are sufficient 
to recover the normal Harman FOM. 



B. Isotropic S, Anisotropic a and ac 

The next case is when the thermopower is taken to be 
isotropic in all directions, S x = S y = S z = S, while the 
electrical and thermal conductivity are still anisotropic. 
Although this situation is unlikely to occur exactly in 
practice, it is theoretically possible. The simplification 
when the thermopower is isotropic is that the Bridg- 
man effect does not occur. Hence, the right-hand side of 
equation (116p is zero, and the temperature distribution 
becomes completely de-coupled from the voltage. From 
equations (|A8[) and ()B2|) . the temperature distributions 
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in the sample are given by 



rpO 



(n ; -n)/c/z 1 



K z ab Vc 2/ 
_ (n, H) J y n mn (a; ; Vi z - K mn ), (25) 



(n t - n) j 



(26) 



/////////// 
t ^ ^^^^^^^^ 



The forms of these expres- 



with if 2 m = *4m 2 + fc 2 2 f 
sions are identical to the corresponding equations for the 
ac voltage, (|24[) and (f2"5| . The key difference is that the 
scale over which the non-uniform part of the tempera- 
ture distribution varies in the z-direction is controlled 
by the thermal anisotropics k xz and k yz , rather than the 
electrical anisotropies. Hence although the mathemati- 
cal forms of the ac voltage and temperature distributions 
may be qualitatively similar, electricity and (thermally 
conducted) heat can flow in very different ways if the 
electrical and thermal anisotropies are different. 

To illustrate this, figure [DJ a) shows the flow pattern for 
the electrical current in a section through a simple sample 
geometry. Figure [2jb) shows the part of the heat current 
due to thermal conduction, jq — Tlj . The sample is ther- 
mally isotropic (k xz = k yz — k — 1), while the electrical 
anisotropy is r xz = r yz = r = 10. While the heat flow 
is confined to a region in the middle of the sample, the 
large electrical anisotropy causes the electrical current to 
spread-out rapidly, leading to an almost uniform flow in 
the bulk. The non-uniform part of the flow occurs only 
in regions with a width of order a/r next to the contact 
faces. Note that the size of the contact and aspect ra- 
tio c/a of the sample are also important here. Smaller 
contacts and aspect ratio make the flows more confined, 
while larger contacts and aspect ratio leads to uniform 
flow over a greater part of the sample. 

Isotropic thermopower also guarantees the absence of 
thermoelectric eddy currents in the sample. This means 
that the dc voltage and temperature distributions in the 
sample are related by the usual expression, \ — X&c + ST. 

Let us explicitly construct the Harman FOM in this 
case. From equation or (f2"5|) (depending on the 
current contact locations), the ac potential difference 
measured between two points (x,y,z) and (x',y',z') is 
Kc = Xac(x',y',z') - Xac(x,y,z) = IRf, where Rf is 
the measured electrical resistance to alternating current 
between the voltage contacts. From equation (|26|) or (125]) 
the temperature difference between the same two points 
is just AT = T{x', y>, z') - T(x, y, z) = -(11,- U)I/K™, 
where Kf 1 is the measured thermal conductance between 
the voltage contacts. Since the relation x = Xac + ST 
holds in this case, the measured dc potential difference 
(including the contribution of the leads) is given by 
Kic = V ac - {Si - S)AT. Hence, the Harman FOM may 
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(b) Jg - nj 

FIG. 2: Plots showing the flow patterns of the electrical and 
heat currents in a cross-section at y — 0.5 through a particu- 
lar sample. The sample has dimensions a = 1, 6 = 1, c = 0.5, 
with current contacts at (0.5a, 0.56, 0) and (0.5a, 0.5b, c). The 
electrical and thermal anisotropies are r = 10, k = 1. The ar- 
rows indicate the local magnitude and direction of the current 
density, with the magnitude capped at 21 jab (2(11 — II;)7/afe 
for the heat current) for clarity, (a) Electrical current j; (b) 
Heat current due to thermal conduction, — nj; 



be written as 



f(Si - S) 2 



(27) 



As expected, the Harman FOM is not a unique number, 
and depends on the voltage contact locations through R™ 
and Kf 1 . However, the physical meaning of the Harman 
FOM survives in this case, since it can be related to a 
measured electrical resistance and thermal conductance. 
Finally, if the electrical and thermal anisotropy in each 



direction are identical (i.e. r xz 



k X z and r yz 



k yz ), 



then R„ 



K mn . Such a situation might occur in a 



metal where the Weidemann-Franz law holds, and the 
electronic contribution to the thermal conductivity dom- 
inates. It is clear from equations (|2"5]) , ([M]) . ([2"5]) and 
([26)1 that the Harman FOM reduces to z H = f(Si - 
Sz^&z/Kz, independent of the details of the current con- 
tact configuration. Essentially, the effects of the small 
current contact modify the electrical resistance and ther- 
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mal conductance in the same way, resulting in no net 
change to the FOM. 



C. Anisotropic Thermopower and Tetragonal 
Symmetry; x = y ^ z 

We will now consider how anisotropic thermopower 
changes the conclusions of the previous section. It is con- 
venient to restrict attention to a case where the crystal 
structure of the material has a tetragonal symmetry with 
x = y ^ z. This provides a minimal level of anisotropy 
which still captures the effects we wish to describe here. 
It will also form a useful basis for the discussion of the 
experiments of Tamura et al in section IVIIII Tetrago- 
nal symmetry means that the sample is isotropic in two 
directions, but anisotropic in the third. We consider 



here the case where the isotropic directions lie parallel 
to the contact faces. In our geometry, the transport co- 
efficients in the x and y-directions are assumed to be the 
same, but take different values to those in the z-direction. 
Hence, the anisotropy parameters satisfy r xz = r„ 



' yx 



= s, k xz — k yz — k. We can extract the m and 
n dependence from the quantities u mni v mn , w mn and Xj , 
each of which is now proportional to uj m n = \/rh 2 + n 2 . 
It is also convenient to define the following combinations 
of parameters 



= z z 1 k 2 — r 2 s(si — s), j3 = si — 1 — z z 



e = -r 2 (s l - s), 



sie, 



S = ai-1, (28) 
? = 0- siS. 



The solutions for each contact configuration are given by 



J 



X' 
S Z T S 



Iz z 
oA 2 



E 

i=l,2 



(-1) 



' + 1 



mI 



M? 



(29) 



\ S z i J a z ab \c 2/ \ 



- a z S z (Hi 
-a z S z (JLi - 



- IT)/k 2 
I1 z )/k z 



IZz 

a 7 d 2 



i=l,2 ^ 2 



-1) 



2+1 



Mi 



*m? 



{x,y,z\\i), (30) 



^ (r 2 + k 2 + r 2 z z (l - s) 2 + {-iy +1 \J (r 2 + k 2 + r 2 z z (l - s) 2 ) 2 - 4r 2 k 2 ^j . (31) 



r 



We note the following important features of this solu- 
tion. First, the voltage and temperature distributions 
both depend on all of the transport coefficients through 
/ij. Previously, the temperature distribution only de- 
pended on the thermal anisotropy through K mn . Second, 
the voltage distributions cannot be written in the form 
X = Xac + ST. These features occur as a direct result 
of the Bridgman effect and thermoelectric eddy currents, 
which couple the flows of heat and electricity together in 
a non-trivial way. 

We will now construct the Harman FOM for current 
contacts on the same side of the sample. The ac potential 
difference again follows from equation (|24[) with R rnn — 
ruj mn , and defines a measured resistance to alternating 
current between the voltage contacts. The temperature 
difference between the voltage contacts is obtained from 
equation 



AT 



Iz, 



(— l) i+1 ' 



<?zS z d 2 [i 2 - [i\ 



(32) 

It is at this point that we begin to see the complica- 
tions introduced by the anisotropic thermopower. Unlike 



before, we cannot formally extract a measured thermal 
conductance from this expression. The reason is that 
the Peltier coefficient is not a unique quantity for the 
sample when the thermopower is anisotropic. One might 
be tempted to write this expression as AT = IP™ I/K™, 
which defines both a measured thermal conductance and 
Peltier coefficient. However, there is little value in this 
approach, since the measured dc potential difference, 
given by 



Vdc 



Iz (— ' 

2 _ 2 («*' + P'rf ) E AAmn (M^m« ) 



Vzd 2 f ff 2 Ml - Mi 



(33) 

cannot be trivially related to either of these quantities. 
Hence the most we can say is that the measured Harman 
FOM is given by equation (fl3l) : 

— — j (a' + P'fJ, 2 ) AA mra (/^aw) 



i=t 2 Ml - Mi 



y^ AA m „ (roj mn ) 



1. 



(34) 
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Similarly, for current contacts on opposite faces of the 
sample we have 



T(Si — S Z ) 2 (T Z ( z' — z\ ab ^ (-1) 



- Z ) + ^2 ( Zz E 2 _ 2 + P'rf ) ^mni^mn) ~ ^ Afl rnn{rU mn ) J 

V i=l,2 ^ 2 mn mn / , . 

z h = — " , (35) 



cd 2 

mn 



which again cannot be related in any simple way to a 
measured electrical resistance, thermal conductance and 
Peltier coefficient. 

We conclude that when the thermopower is 
anisotropic, the Harman FOM becomes even more 
complicated. As well as the dependence on contact loca- 
tions, we cannot even interpret it in terms of 'measured' 
parameters, as we could when the thermopower was 
isotropic. The principle reason for this is the breakdown 
of the relation \ — % ac + ST, which makes it impossible 
to separate the thermoelectric part of the voltage from 
the Ohmic part. In the general case, it would be 
almost impossible to extract any useful information from 
the Harman FOM. However, the degrees of freedom 
associated with the contact locations and the sample 



geometry can be used to simplify the formulae. We will 
consider some of these limits in the next section. 



D. Limit of a Very Thick or Thin Sample 

The solution does take on particularly simple forms 
in certain limits, namely those of a thick sample, and a 
thin sample. For simplicity, we consider a sample with 
a square cross-section, and hence a — b. We take as an 
example the case where both voltage contacts are placed 
on the same face as the current contacts, z = z' = 0. 
The FOM for this configuration is given by 



E(— 1)* +1 , , „/ 2 \ Y^/ w - / ~ / ~ ~ n coth fiiLLi mn c 
— 2\ a + P Mi ) /Xlmn — v mn ){cosmx cos ny — cos trx cos ny) 

i=l,2^-Ml mn Wmn 
Z S H = Z Z ■ ; 1. (36) 

E, . . coth no mn c 
\lmn — f„)(cosmx cos ny — cos rax cos ny) 

r^m.n. 



Now consider the arguments of the hyperbolic functions. 
For c/a — ► oo, the argument is large for all m and 
n, and hence we can approximate cotri^ijCJ mn c/a <~ 1, 
coth rui mn c/a <~ 1 for all m, n. In the opposite limit 
c/a — > 0, terms with small m and n can be approximated 
using coth LiiUj mn c/ a ~ a/(uiL0 mn r), coihruo mn c/a ~ 
a/ (rui mn c). This approximation will fail for large enough 
ra and n, but such terms are irrelevant since the sum is 
always dominated by terms with small m and n. 
Taking the limit c/a^> oo, wc find 



z H -> rz, 



E 



fa 



i^^ -Ml V, 



Vim) - 



1 + k/r + z z (si - l) 2 + z z (si - s) 2 r/k _ J 
v /(l + fc/r) 2 + z z (l-.s)2 



In the opposite limit c/a — » 0, we find instead 

S. 



z H ^z z (si-s) = 



• 2 %?<J X 



ft r 



(38) 



In each of the limits, the sums cancel-out from the Har- 
man FOM, leaving the results depending only on the 
anisotropy parameters. Note that the argument used to 
derive these results did not depend on the contact loca- 
tions on the z = face, or on the size of the current 
contacts. All of the information about these factors is 
lost when the sums cancel-out. In the thick sample limit, 
the FOM is a function of the parameters z z , s, s; and 
r/k only. In the thin sample limit, the Harman FOM 
reduces to the FOM in the x-direction. This is expected, 
since in the limit c/a — ► 0, the FOM can only depend 
on the transport coefficients in the x-y-plane, which are 
isotropic. 
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In reality of course, one never has an infinitely thick 
or thin sample. In deciding whether the system is in 
either regime, the magnitudes of r, /ii and \ii also need 
to be considered. The sizes of the quantities ir/iic/a and 
■nrc/a then determine whether or not the sample is in the 
thick/thin limit. 

There are two other voltage contact locations which 
can be analysed within the general framework presented 
above. We have compiled the results for the Harman 
FOM for these and the previous configuration in table Q] 
Case 2 turns out to be identical to case 1 in the thick/thin 
sample limits, but will give different results in an inter- 
mediate regime. Case 3 is extremely surprising, since the 
Harman FOM is either exponentially large or small in the 
limit c/a — » oo. Again, all of these limits are independent 
of the voltage contact locations on the given faces, and 
of the size of the current contacts. 

A similar analysis can be applied to equation (|35p 
for current contacts on opposite faces of the sample. 
Here, there are only two distinct voltage contact loca- 
tions. The results are also presented in table Q] Case 
4 gives the same results as 1 and 3 in these limits. 
More interesting is case 5, where the voltage contacts 
are on opposite sides of the sample. In the thick sam- 
ple limit, the summation terms tend to zero, leaving 
z° H -> z z (si - l) 2 = f(Si - S z ) 2 (j z /k z . Hence in the 
limit that c/a — > oo, the Harman FOM in this configu- 
ration coincides with the FOM in the z-direction, which 
would be measured in the conventional Harman method. 
In the opposite limit of a thin sample, we instead find 
that the Harman FOM coincides with the FOM in the 
x-direction. 

In fact, case 5 is just the small contact version of the 
conventional Harman method described in section IIIII 
Our analysis suggests that the complications introduced 
by having a small current contact (or an imperfect, large 
contact) can be eliminated by making the sample very 
long and thin. The size of the aspect ratio needed to do 
this depends on a number of factors. As we pointed out 
previously, the particular regime the sample is in will also 
depend on the anisotropy through r, fix and /j,2- However, 
for the thick sample limit of case 5 other factors such as 
the contact locations and current contact size also play 
a major role. The reason is the presence of the constant 
terms in equation (|35[) . For example, let us suppose that 
the sample thickness and anisotropy are such that we can 
make the previous c/a — > oo approximation to Af2 mn . 
Now, if the current contact is sufficiently small, and the 
voltage measured close enough to the current contacts, 
then overall the summation terms may still be larger than 
the constant term. In this case, we would find that the 
Harman FOM would instead be given by equation (f37|) . 
For measurements further away from the current con- 
tacts, the summation terms may become small enough 
to recover the c/a — > oo result, zh — » z z {si — l) 2 . 

The analysis above can be summarised as follows. The 
thin sample limit of all five cases will result in a measure- 
ment of the in-plane FOM, T(Si — S x ) 2 a x /n x . The crite- 



ria for being in this limit are nrc/a <C 1 and /i,rc/a <C 1, 
and in all cases the results only weakly depend on the 
contact locations and sizes. In the thick sample limit 
(defined by nrc/a ^> 1 and y^rc/a ^ 1), cases 1, 2 and 
4 depend only weakly on the contact locations and sizes, 
but the Harman FOM is a reasonably complicated func- 
tion of the anisotropies and the intrinsic FOM. Case 3 is 
also only weakly dependent on the contact locations, but 
gives an either exponentially large or small FOM, and so 
should be avoided. 

In principle case 5 is the best configuration, since it 
offers the possibility to measure the usual Harman FOM 
for the z-direction. Roughly speaking, the criteria for 
being in the thick limit in this configuration with small 
contacts (d 2 /a 2 <C 1) are a 3 /cd 2 yi <C 1 and a 3 /cd 2 r <C 1. 
For large current contacts {d 2 / a 2 < 1), we require instead 
that (a — d)/ HiC <C 1 and (a — d) /re < 1. A compromise 
on the size of the contact would be needed, since it is also 
desirable to place the voltage contacts far away from the 
current contacts. The reasons for this will become clear 
in section I Villi An additional note of caution is that 
although a very long, thin sample will help reduce the 
errors associated with a small/imperfect contact, it may 
increase other sources of error in the experiment. For 
example, errors due to heat loss via radiation would be 
increased by having a very long thin sample. 

VII. LIMIT OF A VERY THIN SAMPLE; VAN 
DER PAUW FORMULAE 

In this section we take a brief diversion to consider 
an interesting application of the theory, which is unre- 
lated to the Harman FOM. Namely, we will prove the 
existence of generalised van der Pauw formulae for thin 
samples. As discussed in the introduction, it has been 
shown by van der Pauw [l2| that two four-terminal re- 
sistance measurements made with point contacts on a 
thin sample of arbitrary shape and isotropic resistiv- 
ity obey a universal relationship. Specifically, consider 
the sample shown in figure [3] Then in the notation 
of section [Vj the four-terminal resistance measurements 
Rabcd = Vabcd/Iab and Rbcda = Vbcda/Ibc obey 
the equation 

exp ( - — Rabcd) + exp ( - ^Rbcda) = 1, (39) 

where t is the thickness of the sample, and p is the resis- 
tivity. 

The proof of this result rests on the following. In 
an isotropic system with point sources/sinks of electri- 
cal current (and neglecting thermoelectric effects), the 
electrical potential satisfies Poisson's equation (c/ equa- 
tion (jT5f ). It is well-known that the solutions of the 
two-dimensional Laplace equation (which is satisfied by 
the potential everywhere except at the points where the 
current is injected/removed) are invariant under confor- 
mal transformation. Therefore, if the electrical potential 
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Case Configuration 



Limit c/a — > oo 



Limit c/a — > 



mi 



s ^ 1 + fe/r + z z (s; - l) 2 + z z {ai - sfr/k _ 



s _^ 1 + k/r + -z z (s; - l) 2 + z z (si - s) 2 r/k _ 

Zh ~* ^/{i + k/ry + z z {i- s y 



i I 



Trc/ a 



Z H — * 2 Z 



l + k/r + z z (si-l) 2 +z z (si-s) 2 r/k 

zh — * ... y 



^(l + k/r) 2 + z z (l-s) 2 



z h -» z z j-^ (si - s) 



X 



Z° H -> - I)' 



o r / \2 

2h -> Z ^( s i ~ S > 



TABLE I: Harman FOM for five different voltage and current contact configurations. A schematic representation of the contact 
configuration is shown in the second column, where current contacts are indicated as squares, and the voltage contacts as 
circles. The sample is orientated such that z = at the bottom. Note that the actual locations of the contacts on the respective 
faces are arbitrary; the configurations shown are purely for illustration. 
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FIG. 3: Schematic representation of the thin lamella to which 
van der Pauw's formula applies. The point contacts labeled 
A, B C and D are placed at arbitrary locations on the edge 
of the sample. 



satisfies a certain relationship in one sample geometry 



(i.e. the van der Pauw formula), then this relationship 
continues to hold in any other sample geometry which 
can be produced by a conformal transformation. In the 
original proof, van der Pauw used an infinite half-plane 
with contacts along the edge to establish equation ([39]) . 



Here, we will investigate whether analogous results 
hold for the temperature and voltage differences mea- 
sured in a thin sample in the Harman set-up. We will 
choose to let the sample thickness tend to zero along the 
y-direction, i.e. b — > 0. The van der Pauw formulae 
which we will derive involve the measured temperature 
and potential differences across the sample. In order to 
properly include the thermopower of the leads, it is con- 
venient to work with the quantity ip — x — S/T = x~ s i 7 1j 
which we refer to as the 'measured' voltage. 

In the thin sample limit, the measured voltage and 
temperature distributions satisfy the following equations 



J 



r — + — + r {si - s)— + {si - 1)- 



dx 2 dz 2 



dz 2 



-{5{x — xi)S(z — z\) — 8{x — X2)S{z — Z2)), 



(40) 



d 2 i] 



d 2 rj 



d 2 ip d 2 ip sil 



s(Si ~ s) )fe2 -( 1 + z * '^te +r2s 7^2 = —{Kx-x 1 )5(z-z 1 )-S(x~x 2 )5(z-z 2 )), (41) 



dz 2 



dx 2 dz 2 a z b 



r 



which are easily derived from equations IT8|) and (fl~9|) . 
Because the measured voltage and temperature do not 
satisfy Poisson equations, the potential and temperature 
differences measured in a sample will not generally satisfy 
van der Pauw type formulae. However, it turns out that it 
is possible to find linear combinations of tp and 77 which do 
satisfy uncoupled Poisson equations. A convenient way 
of finding these linear combinations uses the solution for 
a rectangular geometry which we have found previously. 
Repeating the steps used in the derivation of (|A6[) , we 
find the solution of (|40|) and (|41|) for a rectangular sample 



Iz z 



o-zabc £f 2 /4 - m? 



Ecosmx cos pz , 
~2 Z ,,2rh2 ( C0S 771X2 C0S P Z2 



p z + fl^ffl 



Defining the following function, 



cos mx\ cospzi 



(42) 



4>i(x,z) 



Iz z ■sr-y cos fax cos pz 



a z abc ' p 2 + ajrh 2 



x ( cos mx2 cos pz 2 — cos mx\ cos pz\ I 



(43) 



ip and r\ may be written as 



1 



M2 - Mi 



{a'+Pul) -{a'+P'» 2 )\ f <f> 
{e + 5nl) -{e + 5n 2 ) 



This may be inverted to give 4>± and (f>2 in terms of i\) and 

7]: 



1 



0'e - a'S \ -(e + Stf) (a' + p'fxf) 



T) 



(44) 

Now, it is readily verified that <pi and 02 satisfy the fol- 
lowing equation 

2 d 2 <j>i d 2 <j>i Iz z , . 

-5{x-x 2 )S(z-z 2 )). (45) 

This can be written as a Poisson equation by rescaling 
the x-coordinate via x = fiiX etc: 



d 2 



d 2 



dx 2 dz 



a z b/j,., 



■ {S(x — xi)S(z — zi) 

-5{x-x 2 )S(z-z 2 )). (46) 



Hence, the functions </>i(i, z) and 4>i(x 7 z) are linear com- 
binations of the measured voltage and temperature which 
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satisfy uncoupled Poisson equations. Note importantly 
that the parameters appearing in the transformation (|44[) 
between fa, fa and ip, 77 depend only on the transport 
coefficients, and not on the geometrical properties of the 
sample. Indeed, if this relation is used to replace fa in 
equation (j45|) , then it is straightforward to show that the 
result is simply a particular linear combination of equa- 
tions (|4T)]) and (|4"Tj) . Hence, we conclude that the trans- 
formation of equations (|40|) and (|4Tj) into (|45|) is indepen- 
dent of the rectangular geometry used in the derivation, 
and hence will hold for a sample of arbitrary shape. In 
addition, since the normal derivatives of \ an d V vanish 
individually at the boundaries of the sample, so to do 
the normal derivatives of fa. Hence the boundary con- 
ditions on fa are the same as those for the temperature 
and voltage. 

The analysis of van der Pauw may now be applied di- 
rectly to the fields fa and fa. This leads to van der Pauw 
formulae for differences in fa measured between points on 
the edge of a thin sample of arbitrary shape: 

exp y——j—Afa ABCD ) 

+ cxp (^l A fa BCDA )^l. (47) 
V z x I J 

In the general case, we have been unable to simplify these 
results any further. In the limit that the thermopower 
becomes isotropic, S x = S z = S 7 these formulae reduce 
to 



exp ( ZEVpfj (Vabcd + {S t S)AT ABCD )) 

+ exp {^/^±{y BCDA + (5/ - S)AT BC da)) = 1, 

(48) 



cxp ( - w^m ABCD ) 

+ ^(-w^w ATbcda )- 1 - (49) 

The combination V + (Si — S)AT is recognised as the ac 
voltage, V ac . Hence, the first of these equations reduces 
to the classic form (l39l) of van der Pauw's formula, where 
the resistivity is replaced by the geometric mean resis- 
tivity, (a x a z )'~ 1 l' 2 . Note that the apparent sign differ- 
ence between (14"51) and (|39l) is due to the fact that V and 
V ac are electrochemical potential differences, rather than 
electrical potential differences. Equation (|49|) clearly rep- 
resents the direct thermal analogue of the classic van 
der Pauw formula, where the combination AT/ (Hi — H)I 
is just the thermal resistance between the measurement 
points. 

Finally, we note that the results of other works 0, Hf| 
on the voltage and current distributions in thin samples 



may also be applied to the thermoelectric problem. One 
simply replaces the electrical potential with the field <f> il 
and then uses the transformation (|44|) to relate the results 
to the actual measured voltage and temperature in the 
sample. 



VIII. APPLICATION TO EXPERIMENT 
A. Results 

Having established all of the basic results, we now con- 
sider an application of the theory to some experiments. 
Recently, Kobayashi et al have developed a 'modified' 
Harman method, and used it to measure the figures of 
merit of several anisotropic materials [lH [l6| . The mod- 
ification to the usual Harman method divides into two 
parts. The first concerns the way in which the ac and 
dc voltages are measured, and does not concern us here. 
The second modification is that small current contacts 
are used, and it is their effect that we wish to analyse 
using the formalism developed in this paper. 

We begin with three cases where the modified Harman 
method seems to work well, i.e the Harman FOM ob- 
tained is physically reasonable, and in one case is known 
to be consistent with values from other measurements. 
In the first experiment, the material is Bai^RhsOig, a 
quasi-one dimensional rhodium oxide. Figure a) shows 
the sample and contact configuration used in the mea- 
surement. The sample is orientated such that the cur- 
rent flows along the 6-axis of the material. The key 
features of this configuration are the small contact face 
(a = b = 0.05mm) and large aspect ratio (c/a = 10). 
Since the current contacts extend almost completely over 
the contact faces, the flows of heat and electricity are es- 
sentially uniform. Hence the Harman FOM is expected 
to coincide with the intrinsic FOM in along the 6-axis. 

The sample configurations in the other two cases are 
very similar, and are shown schematically in figure H^b). 
Crystals of two different layered oxides, Bio^sSro^Oa+d 
and BiPbSrCoO, are used, with the current applied par- 
allel to the a-b plane. Both crystals have very small as- 
pect ratios, c/a =0.02-0.03. The layered structure means 
that the largest anisotropy occurs between in and out-of- 
plane directions, with the in-plane directions relatively 
isotropic. Although we have not considered this partic- 
ular combination of sample configuration and anisotropy 
previously, it is straightforward to understand why the 
Harman FOM coincides with the intrinsic in-plane FOM. 
First, because the sample is very thin, the current con- 
tact extends almost completely over the full thickness of 
the sample. This means that the temperature and volt- 
age distributions do not significantly vary along the out- 
of-plane direction. They do however vary along the in- 
plane directions, due to the current contact being much 
smaller than the width of the sample. This leads to non- 
uniform flows of heat and electricity in the in-plane direc- 
tions. However, the transport properties are only weakly 
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anisotropic between the in-plane directions, and hence 
heat and electricity flow along very similar paths. This 
leads to a situation similar to that described at the end 
of section IVII B, where the effects of the small current 
contact cancel-out between the measured electrical re- 
sistance and thermal conductance. Hence, the Harman 
method also works well in this particular configuration. 

In addition to these examples where the Harman 
method works well, two other configurations have been 
investigated by Tamura et al, where the Harman FOM is 
found to substantially larger then the relevant intrinsic 
FOM. Both experiments were performed on BiPbSrCoO 
single crystals, with the contacts faces parallel to the a-b 
plane. To understand the reason for increased Harman 
FOM, we first need a more detailed description of the 
anisotropy in BiPbSrCoO. 

One of the most extraordinary features of this com- 
pound is the degree of anisotropy in the electrical resis- 
tivity between in-plane and out-of-plane directions. The 
out-of-plane resistivity is around 10 4 times larger than 
the in-plane at room temperature, with the degree of 
anisotropy increasing still further at lower temperatures. 
The in-plane resistivit y o f BiPbSrCoO is a few milcm 
at room temperature [26j , making it a relatively poor 
conductor of electricity even in the in-plane direction. 
The large anisotropy with respect to charge transport 
is not seen in the thermal conductivity, with the crys- 
tal structure making the in-plane thermal conductivity 
larger than the out-of-plane by around 7 at room tem- 
perature [27I ]. A final remarkable property of this com- 
pound is that despite the large electrical anisotropy, the 
in-plane thermopower is only larger than the out-of-plane 
by a factor of 2 at room temperature. At room tempera- 
ture, the dimensionless in-plane FOM z a b — 0.029, while 
the large electrical resistivity gives a very small out-of- 
plane FOM, z c = 9.72 x 10 -6 at the same temperature. 

Since the material anisotropy is of quasi-two- 
dimensional form, we make use of the results of sec- 
tion IVI CI The sample is orientated such that the c- 
axis is aligned in the z-direction, with the a and 6-axes 
lying along x and y respectively. The sample has dimen- 
sions in the x, y and z-directions of W x H x L, with 
W = 2mm, H = 2.1mm and L = 0.45mm. All lengths 
will be measured relative to W , giving a sample with 
a 1=1 b = 1, and c = 0.25. The first contact configuration 
considered is shown schematically in figure [5fa). The 
current contacts with a width of approximately 0.065a 
are located at (0.15a, 0.15a, 0) and (0.85a, 0.85a, c), while 
the voltage contacts are located at (0.25a, 0.25a, 0) and 
(0.75a, 0.75a, c). At 180K, the anisotropy parameters of 
the material are r ~ 140, k ss 2.8 and s 1.9. The out- 
of-plane dimensionless FOM, z z , is 1.35 x 10~ 5 at the 
same temperature. In the following we will neglect the 
thermopower of the current and voltage leads, which are 
much smaller than that of the sample in cither direction. 
From our theory, the ac resistance is this configuration is 
predicted to be R ac = 1.05R z q, while the thermoelectric 
resistance is i?dc — -Rac = 12. IT Si/ K z g. This predicts 




FIG. 4: Schematic representations of the two sample config- 
urations used in Kobayashi et al's experiments. The current 
contacts are shown as squares, while the voltage contacts are 
indicated by circles, (a) Configuration for Bai.2Rh80i6, with 
the current directed along the b-axis (b) Configuration for 
Bio.78Sro.403+,j and (Bio.8Pbo.2)2Sr2Co20 y , with the current 
directed along the a-b plane. 



(a) 




FIG. 5: Schematic representations of the two sample config- 
urations used in Tamura et al's experiments. The coordinate 
system is the same as that used in figure [TJ (a) Current con- 
tacts on opposite sides of the sample, (b) Current contacts 
on the same side of the sample. 



the measured Harman FOM to be approximately 12z z . 

Experimentally, the ac resistance, thermoelectric re- 
sistance and Harman FOM at 180K are measured to be 
1.6i? z o, 57TSl/K z0 and 35.5z z respectively. The reason 
for poor agreement between theory and experiment is 
not completely clear, but may be partly due to the con- 
tact configuration employed in this measurement. First, 
the choice of the constant flux boundary condition is ex- 
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pected to produce the largest errors close to the current 
contacts, which is where the voltage measurements are 
made. Second, as we pointed out in the last section, 
this particular contact configuration makes the Harman 
FOM extremely sensitive to the exact positions and sizes 
of the current contacts. To highlight this, we plot in fig- 
ure [5] the FOM measured between points (a — x,a — y, c) 
and (x, y, c), for x and y lying on the line indicated in the 
figure. The maximum value of the FOM occurs between 
the centres of the current contacts, and is around 60z z . 
It is clear that the FOM varies extremely rapidly in the 
vicinity of the voltage contact. The finite size of the volt- 
age contact will also mean that the measured FOM will 
be some average of the region beneath the contact, which 
has not been accounted for here. 




FIG. 6: Harman FOM measured between points (x,y,0) and 
(a — x,b — y,c), for x and y lying on the dotted line indicated 
in the inset. In the experiment the centers of the voltage 
contacts lie at (0.25a, 0.25a, 0) and (0.75a, 0.75a, c). 



The second configuration used in the experiment is 
shown in figure [UJb) . As we pointed out in section IYlD| 
in this configuration the Harman FOM is only weakly de- 
pendent on the contact sizes and locations on the contact 
faces. In fact, the sample turns out to be approximately 
in the thick sample limit, since w/iic/a ps irrc/a ps 110 
and TT[i2c/a ps irkc/a ps 2.2. We can therefore use the 
result for case 3 in table HI which reduces to 

z s H « Zz s 2r -e v{r - k)c/a . (50) 

K 

The exponential factor means that the Harman FOM 
is expected to be enormous, with Zh ~ 1.2 x 10 49 z z ! 
This result occurs because the large electrical anisotropy 
causes the electrical current to flow in a very thin layer 
(with width of order a/r) near to the contact face at 
z = [28| . The ac potential difference measured be- 
tween points on the bottom of the sample at z = c is 
therefore exponentially suppressed, with the theory pre- 
dicting that i? ac = 2 x 10~ 49 i?jo- In contrast, the heat 
current penetrates the sample over a distance of the or- 
der of a/k, and so the temperature difference between the 
voltage contacts is relatively large, with the thermoelec- 
tric resistance predicted to be Rd c — -Rac = 2A8TS Z / 'K z q. 



Experimentally the Harman FOM in this configuration 
is extremely large, although not to the degree predicted 
by the theory. At 180K, the FOM is measured to be 
around 5.5, or equivalently 4 x 10 5 z 2 , which is still an 
extremely large enhancement of the intrinsic FOM. The 
measured ac resistance is 1.1 x 10~ 4 i? 2 o, although the 
origin of this resistance is unclear given that the theory 
predicts that the ac voltage should be effectively zero 
in this configuration. Also puzzling is the magnitude 
of the measured thermoelectric resistance, i?dc — R&c = 
44.2TS z /K z q, which is comparable to that measured in 
the previous configuration. 

In summary, we have used the results derived in section 
IVI CI to qualitatively understand the enhanced Harman 
FOM measured by Tamura et al in BiPbSrCoO. In both 
configurations, we find that the theory underestimates 
the sizes of the ac and thermoelectric resistances in the 
sample. One possible interpretation of the results is that 
the sample is behaving as if the effective electrical and 
thermal anisotropies were considerably smaller than the 
values expected from the resistivities and thermal con- 
ductivities. Alternatively, the uncertainties associated 
with the current and voltage contact locations could also 
explain some of these results. 



B. Physical Origin Of the Enhanced Figure of 
Merit 

In this section we will attempt to understand the phys- 
ical origin of the enhanced FOM predicted by the theory 
in first sample configuration investigated by Tamura et al. 
There are two distinct processes at work. The first con- 
cerns the ways in which heat and electricity flow through 
the sample, while the second involves the anisotropic 
thermoelectric effects. Let us assume temporarily that 
the thermopower is isotropic. In the first contact con- 
figuration the ac voltage and temperature distributions 
are given by equations (f23|) and (|2"5|) respectively, with 
Rmn = ru> mn and K mn = kui mn . Evaluating the ac po- 
tential difference between the voltage contacts gives an 
ac resistance of 1.05i? z oj while the measured thermal re- 
sistance between the same two points is 4.28K~q. Hence 
the measured thermal resistance is increased by a much 
larger factor than the electrical resistance. The reason 
for this is the large difference in the electrical and ther- 
mal anisotropy, which causes electricity and heat to flow 
along very different paths through the sample (cf figure 

From equation (|2"T)) (and taking S = S z ), the Harman 
FOM is evaluated to be 3.8z z . Given that this enhance- 
ment is less than the factor of 12 from evaluating equa- 
tion (|35p . the remaining enhancement must come from 
the effects of thermoelectric anisotropy. To understand 
these effects, we will make some approximations to the 
formulae derived in section IVI CI The large electrical 
anisotropy and small intrinsic FOM suggest that we take 
the combined limit r> 1 and z z <l. To leading order, 
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the quantities fi\ and \i2 are given by 

M 2 «r 2 (l + z 2 (l- S ) 2 ), 
t 4^k 2 (l-z z (l-s) 2 ). 

Substituting these expressions into equations (|30p . and 
making the r » 1 and z z « 1 approximations, we find 
for current contacts on opposite sides of the sample: 



a 2 S 2 (H - U z 



x°(*.y,*) = -^( £ 4)( 

a z ab \c IJ \ 

^((i- ^^ (n,-n z -n z )) ]T a»„(*, v, z\ r) 

- Sx— — — — ^mn(x, y, 2; k) ) , (51) 

mn ' 



r°(i >w ,z) 



/(n ( -n z )c/z i 



n 7 ab 



— 2(( n z - n x )^fi mn (a;,y,2;r) 



(n, -n x )^Q m „(a;,y,z;fc)). (52) 



K 2 d 2 



It is instructive to compare these expressions with those 
for isotropic thermopower, equations (|2"5|) and (|2"B")) . The 
temperature distribution now has an additional term, a 
direct result of the Bridgman effect which allows the cur- 
rent distribution within the the sample to alter the tem- 
perature distribution. This term is very small in this 
case since it is proportional to 1/r. The third term is 
also present in equation (|25p. but in this limit it is the 
in-plane Peltier coefficient which appears in the prefac- 
tor. This is quite surprising, given that because the cur- 
rent contacts are placed on the faces at z = 0, c, and 
the current enters the sample along the z-direction! The 
large electrical anisotropy combined with the Bridgman 
effect are responsible for this. Let us (trivially) rewrite 
the second the term in T° as 

/ 

(EE; - U x ) ^ ^rnn(x, V, Z] k) = 
mn 

t 

((n, - n 2 ) + (n 2 - u x )) *Wz, y, z\ k), (53) 

which we interpret in the following way. The current 
enters the sample in the z-direction, causing a Peltier 
effect proportional to (Hi — II Z ), which is captured by the 
first term. Due to the large electrical resistivity in the z- 
direction, the electrical current turns to flow more easily 
in the x-direction, as shown in figure[2fa). This generates 
a Bridgman effect proportional to (LT 2 — n x ), which we 
associate with the second term. In the limit that r — > 
oo, this process occurs right at the current contact, and 



therefore it is as if the current flows directly from the 
leads into the x-direction. The reverse process happens 
at the other contact. Hence the relative Peltier coefficient 
is (II; — ILc). Note that prefactors of the linear terms still 
involve the out-of-plane thermoelectric coefficients, since 
the uniform flows of heat and electricity resulting from 
them are unaffected by the Bridgman effect. 

The voltage distribution is also modified by the ther- 
moelectric anisotropy. The term proportional to (S z — S x ) 
is interpreted as the Ohmic contribution from the ther- 
moelectric eddy currents in the sample. Again, this term 
is proportional to 1/r and is small in this case. The most 
important feature of the voltage distribution is that the 
final thermoelectric term is proportional to the in-plane 
thermopower. 

The net effect of this is that because the in-plane ther- 
moelectric coefficients appear in the dominant term of 
equation (|5Tj) , the FOM has an enhancement due to ther- 
moelectric effects of roughly s 2 = 3.6. This approxi- 
mately accounts for the difference between the figures of 
merit evaluated from equations ([33]) and (f2"T|) . The actual 
enhancement is less than this, because the the sample is 
not strictly in the r — > oo limit. 



IX. EFFECTIVE FIGURE OF MERIT 

The results of the previous sections have indicated that 
the effects of the small current contact and anisotropy can 
lead to the Harman FOM being larger than the intrin- 
sic FOM of the material in the direction perpendicular 
to the contact face. In the second contact configuration 
investigated by Tamura et al, the measured increase in 
the FOM was extremely large, with the theory predict- 
ing an exponential increase in this case. However, the 
enhancement only occurs for a certain set of voltage con- 
tact locations; there are others for which the Harman 
FOM is actually reduced (see figure [5]). Hence, it is not 
at all clear whether the increased Harman FOM indicates 
an increase in the thermoelectric performance of this sys- 
tem, or is simply an artefact of the measurement. It is 
this question that we shall attempt to address here. 

The central problem with the Harman FOM as defined 
is that it is a local quantity, depending on voltages mea- 
sured between particular points in the sample. In order 
to decide whether the thermoelectric performance is im- 
proved, we need to define a FOM which characterises 
the properties of the whole sample. This can be achieved 
through a consideration of the rate of irreversible entropy 
production in the system. The concept of entropy pro- 
duction as applied to thermoelectric devices is not a new 
idea. Sherman et al [l[ were the first to point-out that 
maximising the efficiency of a thermoelectric device is 
'equivalent' to minimising the rate of unwanted entropy 
production. Specifically, they considered a generator 
made from a thermoelectric material with temperature- 
independent transport coefficients connected to an exter- 
nal load. It was demonstrated that the load resistance 
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which maximises the efficiency also minimises the rate of 
entropy production in the generator divided by the rate 
of entropy production in the load. More recently, Sny- 
der [29[ has introduced the thermoelectric compatibility 
as the quantity which maximises the local efficiency of a 
generator. This quantity may be derived by considering 
the local rates on entropy production. 

Within the framework of irreversible thermodynamics 
[30| . the rate of entropy production is written in terms 
of pairs of conjugate forces and fluxes. For the thermo- 
electric problem, we choose the fluxes to be the electrical 
and heat currents j and jq, for which the appropriate (lin- 
earised) conjugate forces are (1/T)V% and — (1/T 2 )VT. 
The rate of entropy production per unit volume is given 
by 



s = J i dx _ i$i dT 

f dx t f 2 dx t 



(54) 



Using equations (|2|) and j3|) to eliminate jQi and dx/dxi 
from these expressions leads to 



a_1 . . K ik dT 8T 

f PikJi3k+ f2 dxk ; 



(55) 



where the interference term proportional to jkdT/dxi 
vanishes due to the Kelvin/Onsager relation Iljfc = TSki- 
The first term in this expression is the entropy produc- 
tion due to Joule heating, while the second is the entropy 
production associated with heat flow. The total rate of 
entropy production is then the integral of this expression 
over the volume: 



dVS = Sj + S T , 



(56) 



where Sj and St are the total rates of entropy produc- 
tion due to Joule heating and heat flow respectively. 

To gain some intuition about E, let us evaluate it in 
the conventional Harman method. Using equations ([7]) 
and ([5]), we find the total rate of entropy production is 
given by 



S — Ej St = 
I 



I 2 R z0 , (S t ~S z ) 2 I 2 



T 



U ac - (Si ~ S Z )AT) 



K z0 
T 



(57) 



Hence the total rate of entropy production is simply re- 
lated to the power obtained from the measured (two ter- 
minal) voltage and the current. The key observation from 
this analysis is that the FOM of the sample may be writ- 
ten in terms of St and Sj: 



St = T(S l - S z f 
Sj R z qK z0 



(58) 



The advantage of this expression compared to the Har- 
man definition (|13|) is that St and Sj are quantities 
which take single values for a sample with a particular 



configuration of current contacts. The freedom with re- 
spect to the locations of the voltage contacts has been 
eliminated. 

Let us now apply this result in the regime of small cur- 
rent contacts and isotropic thermopower. In appendix 
[Cl we show that equation (|57|) continues to hold, but 
where the potential and temperature differences are re- 
placed by contact averaged values. The contact aver- 
aged potential/temperature difference is defined as the 
average potential/temperature difference between points 
lying on the current contacts. The rate of entropy pro- 
duction due to Joule heating is related to the contact 
averaged ac voltage via TEj = I(V ac ). We can there- 
fore use this relation to define an effective electrical re- 
sistance for the sample, R e g = (V ac }/I = TT,j/I 2 . This 
is the resistance which is 'seen' by the external voltage 
source which drives the current through the sample. Sim- 
ilarly, the rate of entropy production due to heat flow 
is related to the contact averaged temperature differ- 
ence across the sample via TEt = — (Si — S)I(AT). 
Hence an effective thermal conductance can be defined 
as K cS = -{Ui - n)//(AT) = T 2 E T /(AT) 2 . Finally, 
the effective FOM is given by 



St = T(Si - S) 2 
Sj R e gK c f[ 



The effective FOM is related to an effective electrical re- 
sistance and thermal conductance, which reflect the flows 
of heat and electricity through the whole sample. 

The true value of defining the FOM in terms of the en- 
tropy production in the Harman set-up is only seen when 
the thermopower becomes anisotropic. For anisotropic 
thermopower, the individual rates of entropy production 
due to Joule heating and heat flow are not formally re- 
lated to the contact averaged ac potential difference and 
temperature difference respectively. However, the rela- 
tion between the total rate of entropy production and the 
contact averaged dc potential difference, TE = J(Vdc), 
does still hold. The problem previously was that we did 
not know how to divide the dc voltage into an Ohmic 
part and a thermoelectric part; the rate of entropy pro- 
duction shows how this division should be made. If we 
define (V dc ) = IR cS + S cS AT, and AT = TS cS /K cS , 
then the effective resistance, thermopower, thermal con- 
ductance and FOM are defined by 



ReS = 



TS, 
I 2 



Sr.fi = 



ZcS 



St 

s7 



TS T 
/(AT) ' 
= f(S eS ) 2 

ReSK c S 



K eS = 



(AT) 5 



where Sj and St are calculated from equations (|29|) or 

dm. 

As an example of using these formulae, we will cal- 
culate the effective transport quantities for the samples 
used in Tamura et al's experiments. In the first con- 
figuration (contacts on opposite sides), we find that the 
effective resistance is 1A3R z q, close to R z q as expected on 
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account of the large electrical anisotropy. The effective 
thermal conductance is 0.044AT z0 , which is substantially 
reduced on account of the weak thermal anisotropy and 
small current contact. The effective thermal resistance 
is therefore increased much more than the effective elec- 
trical resistance. The effective thcrmopower is 1.87 S z , 
which is very close to S x as predicted from the Harman 
FOM. Finally, the effective FOM is 55.3z z , which is con- 
sistent with an average value over the current contact 
(see figure [S]). 

In the second configuration, we find the effective re- 
sistance, thermopower and thermal conductance to be 
0.42i? z0 , 1.895* and 0M7K zO . It is more difficult to in- 
terpret the results from this configuration, since the 'nor- 
mal' results for large contacts cannot be obtained in this 
case. The most important feature is that the effective re- 
sistance is not exponentially small, unlike the measured 
resistance on the top of the sample. This is testament 
to the fact that the effective resistance corresponds to 
the physical resistance of the sample. It is however sub- 
stantially smaller than i? z o, simply because the electri- 
cal current is not forced to flow in the z-direction where 
the resistivity is largest. The effective thermal conduc- 
tivity is still substantially smaller than K z q, while the 
effective thermopower is again very close to S x . The ef- 
fective FOM is much larger in this configuration, with 
z e s = 178z z , primarily due to the lower electrical resis- 
tance. 

We conclude that the effective electrical resistance and 
thermal conductance describe the ways in which electric- 
ity and heat flow through the sample, while the effec- 
tive thermopower captures the thermoelectric processes. 
All three quantities include the effects of eddy currents 
and the Bridgman effect in an appropriate way. How- 
ever, we have not yet shown that the effective FOM is 
in any way related to the efficiency of a thermoelectric 
device. Suppose we used the sample configuration shown 
in figure H^a) to create a simple thermoelectric genera- 
tor. Hot and cold reservoirs would be connected over 
the current contacts, generating a unform thermoelec- 
tric voltage between the current contacts which could 
then be applied to an external load. Our analysis of the 
Harman method implies that due to the anisotropy the 
thermal resistance will be increased much more than the 
electrical resistance. Further, the thermoelectric voltage 
will be increased due to the eddy currents and Bridg- 
man effect. We would therefore expect the thermoelectric 
performance to be improved compared to the situation 
where the reservoirs are connected across the full area of 
the sample. Since the effective FOM of the sample in this 
configuration is increased relative to z z , we conclude that 
the effective FOM is at least qualitatively related to the 
thermoelectric performance of the sample when viewed 
as a thermoelectric device. 

Although we have shown the effective FOM is en- 
hanced compared to the intrinsic FOM in one direction, 
we have not yet compared it with the intrinsic figures 
of merit in the other coordinate directions. For exam- 



ple, in the first configuration of Tamura et al's experi- 
ment, we found that effective FOM was increased by 55.3 
compared to z z . However, the intrinsic FOM in the x- 
direction for this material is z x — (r 2 s 2 /k 2 )z z — 9025z z ! 
Therefore, although the effective FOM is substantially 
enhanced compared to the intrinsic FOM in one direc- 
tion, it is still much smaller than the intrinsic FOM in 
the other. Recall that if the current is passed into the 
sample in a particular direction through large contacts 
(and with current leads of zero thermopower), then the 
effective FOM is just the intrinsic FOM for that partic- 
ular direction. Therefore, one can only really claim that 
the thermoelectric performance is improved with small 
contacts if the effective FOM exceeds all of the intrinsic 
figures of merit. Otherwise, it would always be better to 
pass the current though large contacts in the direction 
with the highest intrinsic FOM. We will explore this is- 
sue in the next section, by establishing some bounds on 
the effective FOM. 



X. BOUNDS ON THE EFFECTIVE FIGURE OF 
MERIT 

The aim of this section is to establish some bounds 
on the effective FOM. It is convenient for this analysis 
to set the lead thermopower Si equal to zero. None of 
the results we present rely on this assumption, but the 
analysis becomes considerably simpler with it [3l|. We 
will prove that in an anisotropic system with any con- 
tact configuration, the effective FOM is always bounded 
from above by the largest intrinsic FOM of the material. 
Here, we will prove the result in the case of tetragonal 
symmetry and isotropic thermopower in order to illus- 
trate the method. This argument is then extended to 
include anisotropic thermopower in appendix [D] The 
simplest way to prove this result is to return to equa- 
tions (|2"Tj) and (j2"2"|) . Restricting to tetragonal symmetry 
and isotropic thermopower, and also setting Si — 0, the 
ac voltage and temperature distributions may be written 
as 



A, 



^ cos rhx cos hy cos pz, (59) 



mnp mn c 

T= — — — mnp — - cos rhx cos ny cos pz . (60) 
k z ^ k z LO^ n + p z 

mnp " tlL 

The quantities A mnp contain all of the information about 
the boundary conditions associated with the current con- 
tacts. For example, for current contacts on the face at 
z = 0, we would have A mnp = (j mn - v mn )/cd 2 . The 
proof does not require any assumptions about the posi- 
tions or sizes of the current contacts, with all of the rel- 
evant information subsumed into A mnp . The important 
point is that the information about the current contacts 
enters both Xac and T in the same way through A mnp . 

Using these expressions the rates of entropy production 
Sj and T,t may be derived. The effective FOM is found 
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to be 



i 



A 2 

mnp 



^ k 2 uj 2 



f>2 



ZoS 



E 

mnp 



A 2 

mnp 



(61) 



T 2 UJ 2 -\- V 2 



We now compare corresponding terms (i.e. those with 
the same values of m, n, and p) in the numerator and 
denominator. Note that every term in the numerator and 
denominator of this expression is positive. For r > k, it 
is clear that every term in the numerator is greater than 
or equal to the corresponding term in the denominator. 
Equality holds for u> mn = 0. Therefore, the sum in the 
numerator is greater than or equal to the sum in the 
denominator. Hence for r > k it follows that z e g > 
z z . Conversely, for r < k every term in the numerator 
is less than or equal to the corresponding term in the 
denominator, and z e g < z z . 

The next step is to investigate the relative sizes of 
z e f£ and z x . Noting that for isotropic thermopower 
z x = (i" 2 /k 2 )z Zl equation (|6 1 1) can be rewritten as 



E 



A 2 

mnp 



p 2 /k 2 



Z c ff 



E 



A 2 

mnp 



(62) 



u mn +P 2 /r 2 



The analysis above can now be repeated for z x . For 
r > k every term in the numerator is less than or equal to 
the corresponding term in the denominator, and conse- 
quently z c ff < z x . Conversely, for r < k every term in the 
numerator is greater than or equal to the corresponding 
term in the denominator, and so z e g > z x . 

A summary of the various cases is given below: 

for r/k > 1, z z < z cB < z x ] 
for r/k < 1, z x < z cS < z z ; 
for r/k = 1, z z = z cS = z x ; 

where the final case trivially follows from equations (|6"Tj) 
and ([62]) . We conclude that when the thermopower is 
isotropic the effective FOM is always less than or equal to 
the largest intrinsic FOM, in any contact configuration. 

The extension of this proof to anisotropic thermopower 
may be found in appendix [Dj The analysis is very sim- 
ilar to that presented here, but with several additional 
complications. For anisotropic thermopower the intrinsic 
figures of merit are related by z x = (r 2 s 2 /k 2 )z z . The pa- 
rameter which controls the relative sizes of z e g, z z and z x 
is r\s\/k. The modulus signs are needed since s — S x /S z 
could be either positive or negative. The various cases 
are: 



for r\s\/k > 1, z z < z e g < z x or z e g 
for r\s\/k < 1, z x < z e g < z z or z e g v z x v z 2 ; 
for r\s\/k = 1, z e g < z z and z e g < z x . 



< z z < z a 

< z x < z z 



Again, it is clear that z e g is always less than or equal 
to the largest intrinsic FOM. The added complication 
when the thermopower is anisotropic is that it is much 
harder to define the regime where each intrinsic FOM is 
larger than z c ff. This regime partly depends on other 
factors such as the current contact locations and sizes. It 
is also possible for z e g to be simultaneously smaller than 
both intrinsic figures of merit, which was not possible for 
isotropic thermopower. 

In summary then, we have shown in this section that 
for a sample with tetragonal symmetry and an arbitrary 
arrangement of current contacts the effective FOM is al- 
ways bounded from above by the largest intrinsic FOM 
of the material. Therefore, if the maximum FOM is the 
only consideration, it is always better to pass the current 
through large contacts in the direction where the intrin- 
sic FOM is highest. It should be possible to extend this 
analysis beyond tetragonal symmetry to the most gen- 
eral case of anisotropy. We also note that our results are 
somewhat similar to those of Bergman et al [32] ]. who 
proved that for a composite thermoelectric material the 
overall FOM cannot exceed any of the individual figures 
of its constituents. 



XI. CONCLUSION 

We have investigated the effects of small current con- 
tacts and anisotropy on the measurement of the FOM 
using the Harman method. Our model captures the non- 
trivial effects associated with thermoelectric anisotropy 
in the form of thermoelectric eddy currents and the 
Bridgman effect. This introduces a self-consistency be- 
tween the temperature and voltage distributions. Using 
our solutions, we investigated what would be measured 
in the Harman method with small current contacts. We 
found that the Harman FOM, zh 1 is generally very dif- 
ferent from the intrinsic FOM of the sample, except in 
a few special cases. When the thermopower is isotropic, 
the Harman FOM may still be interpreted in terms of a 
measured electrical resistance and thermal conductance. 
This interpretation is no-longer possible when the ther- 
mopower is anisotropic, primarily due to the breakdown 
of the relationship between the dc, ac and thermoelec- 
tric voltages. We also investigated how geometrical fac- 
tors can also change the Harman FOM. The small con- 
tact version of the conventional Harman method (cur- 
rent and voltage contacts on opposite sides of the sam- 
ple) is particularly sensitive to the anisotropy, geome- 
try and current contact sizes and locations. We believe 
that this is the main reason reason for the poor quantita- 
tive agreement between the theory and the Harman FOM 
measured in the first configuration of Tamura et al's ex- 
periments, although we cannot rule-out some other non- 
trivial effects not captured by this theory. A key test of 
the theory would be to measure the Harman FOM in one 
of configurations 1 , 2 or 4 of table U in the thick sample 
limit. The Harman FOM in these cases depends on all of 
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the anisotropy parameters, but is only weakly dependent 
on the contact sizes and locations. 

We also extended the analysis of van der Pauw for 
thermoelectric measurements thin samples of arbitrary 
shape. We showed that there exist van der Pauw-type 
formulae relating particular linear combinations of the 
potential and temperature differences measured between 
points on the edge of the sample. 

The work on thermoelectric measurements in 
anisotropic systems is not yet complete. Further analysis 
is required to determine the effects of heat losses, which 
were neglected in our calculations. Overall though, the 
main message from this part of the work is that great 
care should be taken in interpreting the FOM measured 
using the Harman method, particularly in anisotropic 
systems. 

In the second part of the paper we explored whether 
the increased Harman FOM measured in certain sample 
configurations actually indicates an improvement in the 
thermoelectric performance. We argued that an effective 
FOM could be defined from the rates of entropy pro- 
duction occurring in the sample in the Harman set-up, 
z e g = St/Sj. Unlike the Harman FOM, this quantity 
takes a single value for a sample with a given configu- 
ration of current contacts. We identified two effects of 
small current contacts and anisotropy which can cause 
the effective FOM to be enhanced. First, the flow of 
heat can confined relative to the flow of electricity, lead- 
ing to the effective thermal resistance being increased 
more than the effective electrical resistance. Second, the 
effective thermoelectric power can also be enhanced by 
the Bridgman effect combined with thermoelectric eddy 
currents. However, we also proved that for a sample with 
small current contacts and whose transport coefficients 
have tetragonal symmetry, the effective FOM is always 
bounded from above by the largest intrinsic FOM of the 
material. 

Despite this fact, the use of anisotropy to improve ther- 
moelectric performance in the manner indicated in this 
work may still be valid. It is important to remember 
that the intrinsic FOM only determines the efficiency of 
a thermoelectric device under certain assumptions (as de- 
scribed in the introduction), and is generally used as a 
guide for selecting materials as potential thermoelectrics. 
Other factors such as the mechanical properties of ther- 
moelectric materials or manufacturing issues are also im- 
portant considerations in developing thermoelectric de- 
vices. The primary reason for the enhanced effective 
FOM in an anisotropic system with small current con- 
tacts is that heat and electricity can be made to flow 
along different paths. This ability to independently con- 
trol the flows of heat and electricity in anisotropic sys- 
tems is one of the main findings of this work, and may 
itself have useful applications. 
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APPENDIX A: SOLUTION FOR CURRENT 
CONTACTS ON THE SAME FACE OF THE 
SAMPLE 

In this configuration, both contacts are located on the 
face at z — 0, and hence we set z\ = = z-i in equations 
(f2~Tj) and (|22[) . At this stage, we can also remove the 
restriction to point contacts. To do this, we 'broaden' 
the delta-functions in the x and y-directions into top-hat 
functions of width d; 5(x — xi) — > (l/d)9(d/2— \x — Xi\) 
etc. The point contacts results are recovered in the limit 
d — > 0. The Fourier representation of the contact located 
at (xi,yi,zi) is given by 



f 

rf2 



Q(d/2 - \x - Xl \)Q(d/2 - \y - Vl \) 



-y 

d? ^ 



v mn cos rhx cos fiy; (Al) 



v oo — —r- 

ab 

d 2 mird 



b (rmr) 2a 
d 2 rmd 



a (mr) 2b 
4 



sm — — cos nyi , 



mud . nnd 
sm — — sm — cos mx \ cos nyi . 



(m7r)(n7r) 2a 2b 



(A2) 



The second contact, located at {xi,yi,z-i) has a similar 
Fourier representation with coefficients T mn , which may 
be obtained from v mn with the replacement 1 — > 2. 

A considerable simplification of (j2Tj) and (f22|) can be 
obtained by performing the sums over p. To do this, we 
define the following anisotropy parameters 



(A3) 
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\a z t>t n z ' 



"7,,,, : : ( — -5r H z z - jm 

( r xz s xz "I" k xz z z )m + {fy Z Sy Z + k yz z z )h 



(A5) 

We have also defined the intrinsic FOM in the z-direction, 
z z = TS z a z / k z . In terms of these parameters, x s an d 
T s may be written as 
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map 
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(A7) 



where we have also defined a dimensionless thermopower of the current leads, si — Si/S z . In going from the first to 
the second line, the denominator has been factorised into the form (p 2 + X 2 )(p 2 + A 2 ,), and then the expression split 
into partial fractions. This allows the sum over p to be performed, giving the final result 



X 3 \ _ 1 Z * W ( ™ln - SlV 2 mn + (81 ~1-Z- Z ')X 2 \ 

S,T> 2^1. A 2 -A 2 [ v 2 mn - Sl u 2 nn + ( Sl -l)X 2 JA ro „( W A), (A8) 

ran i— 1,2 ^ ■ 



where 

^■mn{*C-> ?/, Z 1 Aj) 

. coshAi(z — c) _ _ . . 

(7m»-M-; — r-n cosmxcosny. (A9) 

A.; Sinn XiC 

Note that in the original representation (i.e. before sum- 
ming over p) all terms with m = n = vanished indepen- 
dently of p. Hence, the term with m = n = is excluded 
from the sum in equation (j A8|) . 

Note that in having performed the sum over p, the 
boundary conditions satisfied by the functions \ an( i V 
have changed. The solutions now satisfy homogeneous 
versions of equations (fl8|) and (Tl9|) i.e. without the 
source/sink terms on the right-hand side. The bound- 
ary conditions on the faces at x = 0, a, y = 0, b and z = c 
remain unchanged, but the boundary conditions on the 
face at z — are altered to read 
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These boundary conditions correspond to a constant flux 
of electrical current (and Peltier heat) into the contact 



at (#1,2/1,0) and out of the contact at (x2, 2/2,0). Indeed, 
if we had started with this boundary condition initially, 
along with the homogeneous forms of equations (jT5J) and 
([12]), then equations (|A8|) would have been obtained di- 
rectly. 

APPENDIX B: SOLUTION FOR CURRENT 
CONTACTS ON OPPOSITE FACES OF THE 
SAMPLE 

In this configuration, one contact is placed on the face 
at z = 0, while the other is placed on the face at z = c. 
Hence, we simply set z\ = and zi — c in equations (|21|) 
and (122|) . The term proportional to j mn now acquires a 
prefactor of (— l) p . This leads to an important difference 
with the previous case, namely that not all of the terms 
with m = n = vanish. The terms with m = n = are 
given by 
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Rewriting the remaining terms in the same way as before, 
and performing the summation over p, we find the final 
result for contacts on opposite faces 
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A; sinh A;c 



cos mi cos ny. (B3) 



APPENDIX C: ENTROPY PRODUCTION AND 
CONTACT AVERAGING 

Here we will prove the results quoted in the main text, 
which relate the rates of entropy production to contact 
averaged temperature and potential differences. We will 
prove the results for isotropic thermopower. In addition, 
we will use a sample configuration with the finite-size 
current contacts on the same side of the sample at z = 0. 
This choice is arbitrary, but the analysis is simpler to 
follow with a definite example. Given these restrictions, 
the voltage and temperature distributions in the sample 
satisfy 
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All surface terms vanish due to the insulating boundary 
conditions on the surfaces of the sample. The quantity 
X — ST is recognised as the ac voltage, %ac- Substituting 
for A(x, y, z), this expression may be re- written as 
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The quantity in square brackets is just the average ac 
potential difference between the current contacts. Hence 
the total rate of Joule heating in the sample, TEj is 
related to the average power supplied by the external 
current source, I (Vac). 

Similarly it can be shown that the rate of entropy 
production due to heat flow, Er, is related to the contact 
averaged temperature difference: 

dT\ 2 



mv, /"W ( dT \ 2 ( dT \ 2 ( dT 
= J \ x \~dx) + Ky \dy) + <lTz 



= -I(St-S)(T) 



(C6) 



Hence the total rate of entropy production is related 
to the contact averaged dc potential difference, TE = 
I(V dc ) = /((Vac) - (Si - 5) (AT)). 



A(x,y,z) = ^(9(d/2 -\x- x 1 \)Q(d/2 - \y - yi \) 
- Q(d/2 -\x- x 2 \)G(d/2 -\y- y 2 \))S(z). (C3) 

Consider the rate of entropy production due to Joule 
heating, Ej. Substituting for fa using equation (fT]), inte- 
grating by parts, and using equation (jCljl provides 

TEj = f dv(^ + ^- + ^) 

= - I dV( X -ST)IA(x,y,z) (C4) 



APPENDIX D: BOUNDING THE EFFECTIVE 
FIGURE OF MERIT; ANISOTROPIC 
THERMOPOWER 



We will prove that in the case of tetragonal symmetry 
and anisotropic thermopower the effective FOM is always 
bounded from above by the largest intrinsic FOM of the 
material. Using equations (f2Tj) and (|22|) and setting Si = 
0, the effective FOM is found to be 
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Zcff = Zz 



j2^(kW mn+ p 2 )(r 2 scu 2 mn+ p 2 ) 2 

mnp U rnnp 
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— D 2 

mnp mnp 



mnp "mnp 

r 2 LU 2 mn (k 2 LU 2 mn + (1 + (1 - s)z z )p 2 y +f({k 2 - r 2 z z s{\ - s))u 2 mn +p 2 ^ 



(DI) 



where D mnp — (r 2 uj 2 nn +p 2 )(k 2 uj 2 an +p 2 ) + r 2 z z (l — s) 2 oj 2 nn p 2 . The form of equation (|D1|) indicates that z c g is always 
positive. Expanding the numerator and denominator, this equation has the form 



' A 2 

E mnp / 

mnp rnnp 
Zett = Z z 



' A 2 

4 ~2 i 2 ~4 , ~6\ -"-mnp ti% , l 4 ~2 , / 2 -4 , / ~6\ 

n 2 u mn p + n 3 w m „p + n 4 p ) 2^ Jyi — K w m » + n 2^mnV + "a^mnP + n 4 p ) 

mnp TfiTip 
= ' .12 



'A 2 

mnp mnp 



' A 2 

E^mnp / j 6 i j 4 ~2 , j 2 -4 , j -6 
Jy2 [dl^mn + d 2^mnP + d^mnP + «4P 

™ rnnn 



mnp Tnnp 



We begin by establishing the regimes in which z e g is 
larger, smaller or equal to z z . To do this we consider the 
relative sizes of the coefficients of like powers of w mn and 
p in the numerator and denominator of the first equality 
in equation (|D2p i. e. m and di . The second equality will 
be used later to establish the size of z e ff relative to z x . 
The factors of -A^ / may be ignored since they 
appear in both the numerator and denominator. The 
coefficients di of the terms in the denominator are all 
positive, as are ni and 714. The coefficients 712 and 713 
can be negative for s < 0, which occurs if the intrinsic 
thermopowers have opposite signs. If we can show that 
Tii/di > 1 for all i, then the numerator of the first equal- 
ity of equation (|D2j) is definitely greater than or equal to 
the denominator, and z e ff > z z . Similarly, if we can show 
that rii/di < 1 for all i, then the numerator of equation 
(|D2[) is definitely less than or equal to the denominator, 
and z e g < z z . This is essentially the strategy we will use 
in this proof. 

Taking first 0(uj^ ln ), we have n\/d\ = T 2 , where T 2 — 
r 2 s 2 /k 2 . Hence it is clear that if T 2 > I, rii/di > 1, while 
for r 2 < 1, n\jd\ < 1. The terms of 0(p 6 ) are also easily 
bounded. We find that n^/d^ = 1, independently of the 
value of r 2 . In figure E£a) we indicate the behaviour of 
ni/di and 774/^4 in the T 2 -s plane. The reason for doing 
this will become clear shortly. 

Next we deal with the terms of 0{bJ^ nn p i ): 



"2 
~2 



s 2 + 2T 2 + 2T 2 z z {\ - s) 2 + T 4 z 2 (l - sf 



(D3) 



Unlike the previous two terms, n 2 /d 2 is not simply a func- 
tion of r 2 , but also depends explicitly on s and z z . Not- 
ing that r 2 is necessarily positive, the boundary where 



n 2 /d 2 = 1 is given by 
_ 2 l-,s + z,(l-s) 2 



l-z 2 {l-s) 2 

yj(l-s + z z {l- s) 2 ) 2 + s 2 (l - z 2 2 (I - s) 2 ) 
l-z 2 (l- S ) 2 



(D4) 



In figure EJb) we plot this expression as a function of s 
for a fixed value of z z . The positive part of this function 
is contained within the two asymptotes at z = 1 — l/z z 
and 1 + l/z z . It can be shown that the minimum value 
this function is T 2 — 1 and occurs at s = 1 independently 
of the value of z z . For T 2 lying above this line we have 
n 2 /d 2 > I, while for T 2 lying below the line, n 2 /d 2 < 1. 

The terms of C(a;^ ln p 4 ) are analysed in exactly the 
same way: 



d 3 



s 2 + 2T 2 s 



2s 2 + r 2 + 2r 2 z z (i - S ) 2 + r 2 z 2 (i - s ) s 



The boundary where n^/ds = 1 is given by 



2s - 1 - 2z z (l - s) 2 - z z (l - s) 2 



(D5) 



(D6) 



This function is plotted in figure 0c) . The positive part 
of this function is contained within the two asymptotes 
at z = (1 + z z )/(2 + z z ) and 1 + 1/ z z . As before, it can be 
shown that the minimum value of this function is T 2 = 1 
and occurs at s = 1 independently of the value of z z . For 
r 2 lying above this line we have 773/^3 > 1, while for T 2 
lying below the line, n^/dy, < 1. 

We can now establish the relative sizes of z z and z c ff- 
The most important regime to identify is where z c g is 
definitely less than or equal to z z . The condition for this 
to be true is that rii/di < 1 for all i. An examination 
of figures [7[a)-(c) reveals that this condition is met for 
r 2 < 1, independently of s. There is also region where 
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n\/di > 1 
n'i/di — 1 


„ ij i 

n^/ct^ — 1 

n'Jch < 1 


1 








rii/di < 1 


n^jd\ — 1 




n'l/di — 1 


n' A /d± > 1 



(a) 0«J and 0(p 6 ) 




(b) 0{^ mn f) 




(c) 0(^„p 4 ) 

FIG. 7: Plots indicating the relative sizes of the coefficients 
appearing in the expression for z c g. Solid lines are curves 
where rii/di — 1, while dotted lines indicate asymptotes to 
these curves, (a) Coefficients of terms of C(w^ l „) and 0(p 6 ); 
(b) Coefficients of terms of C^w^p 2 ); (c) Coefficients of 
terms of 0{uJ 2 rln p 4 '). 



rii/di > 1 for all i, meaning z c g is definitely greater than 
or equal to z z . Finally, there is a region where definite 



statements cannot be made about the relative sizes of 
z e ff and z z . In this region, rii/di > 1 for one or more z, 
and rii/di < I for the others. Here, the overall size of z c s 
relative to z z depends on a competition between all of the 
terms in the sum, and will also depend on the quantities 
A m np which encode the current contact positions and 
geometry. 

Having established the bounds on the effective FOM 
relative to z z , we now need to repeat the entire procedure 
for z x . For anisotropic thermopower the relationship be- 
tween z z and z x is z x — (r 2 s 2 /k 2 )z z = T 2 z z . We now use 
the second equality in equation (|D2[) . where n! i = rii/T 2 . 
At ©(w^) we find n' x jd\ = 1 independently of T 2 , while 
at C(p 6 ) we have n' 4 /d4 = 1/T 2 . The various cases for 
n'i/di and n'Jd^ are shown in figureJTJa). For 0(uim n p 2 ) 
the curve on which n' 2 /d2 = 1 is given by 

l + 2z,(l- S ) 2 



r 2 



~ 2^ 2 (l - ^) 2 

+ 2z z (l - s) 2 ) 2 + Asz 2 (2 - a)(l - s) 5 
2^ 2 (1 - s) 2 



, (D7) 



which is sketched in figure Wft>) . This function is only 
positive for < s < 2, and takes its maximum value at 
s = 1, r 2 = 1 independently of the value of z z . Finally, 
for 0(ui 2 nn jfi') the curve on which n'^/d^ = 1 is given by 



I + 2z 2 (l-s) 2 + 2 2 (l-s) 2 
yj(s- s 2 ) 2 + s 2 (l + 2z z (l - s) 2 + z 2 (l - s) 2 ) 



l + 2z z (l-s) 2 + z 2 (l~s) 2 



(D8) 



which is plotted in figure E^c) . This function is positive 
for all s, and also takes its maximum value at s = 1, 
T 2 = 1. 

We can now set bounds on the size of z e g relative to z x . 
Again, the most important region to define is where z e ff 
is definitely less than or equal to z x . The condition for 
this is that n'Jdi < 1 for all i. A consideration of figures 
[2a)- (c) reveals that this condition is met for T 2 > I, 
independently of s. As for z z , there is also a region where 
z e fi is definitely greater than or equal to z x , and another 
where definite statements cannot be made. 

We now have all of the information required to prove 
the result that the effective FOM is bounded from above 
by the largest intrinsic FOM. The crucial statement is 
that the T 2 -s plane may be divided up into two regions. 
For r 2 < 1, z e s is definitely less than or equal to z z , 
independently of s. For T 2 > 1, z e g is definitely less 
than or equal to z x , again independently of s. This is 
sufficient to prove the result. At any point in the T 2 - 
s plane, we can say with certainty that z e ff is less than 
or equal to the largest intrinsic FOM. The fact that the 
relationship between z c g and the other intrinsic FOM 
may be indeterminate is irrelevant. Finally, we note that 
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this proof fully reproduces the results of section [X] for line s = 1 in the T 2 -s plane, 
isotropic thermopower, which follow by restricting to the 
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